


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1972 


Design considerations for an experimental 
heated flat plate stability study. 


Lowell, Robert Leroy. 


Case Western Reserve University 


http://hdl.handle.net/10945/16065 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (q D U DLEY research materials and institutional publications created by the NPS community. 
«ist | Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


w KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


DESIGN CONSIDERATIONS FOR AN EXPERIMENTAL 
HEATED FLAT PLATE STABILITY STUDY 


Robert Leroy Lowell 











DESIGN CONSIDERATIONS FOR AN EXPERIMENTAL 


HEATED FLAT PLATE STABILITY STUDY 


by 


ROBERT LEROY LOWELL, JR. 


Submitted in partial fulfillment of the requirements 


for the Degree of Master of Science 


Department of Fluid, Thermal, and Aerospace Sciences 
School of Engineering 


CASE WESTERN RESERVE UNIVERSITY 


- 


January 1972 


LIBRARY A 
HAVAL POSTC 


RADUATE SCHOON 
MONTEREY, CALIF.. 93940 —— 


CASE WESTERN RESERVE UNIVERSITY 


GRADUATE STUDIES 


We hereby approve the thesis of 


ROBERT LEROY LOWELL, JR. 
candidate for the ES, 


degree. 





DESIGN CONSIDERATIONS FOR AN EXPERIMENTAL 
HEATED FLAT PLATE STABILITY STUDY 


Abstract 


by 


ROBERT LEROY LOWELL, JR. 


The design and construction of an apparatus for experimentally 
examining the stability of a uniform, parallel water flow over a 
heated, constant-temperature, sharp-edged, flat plate to a 
mechanically generated disturbance is discussed. Consideration 
is given to the effects of the streamwise pressure gradient, test 
section turbulence intensity level, uniformity of the flow and 
orientation of the plate in it, material considerations for the 
plate, maximum utilization of the existing facility, heat genera- 
tion required to maintain a constant temperature plate for different 
flow speeds, and production of a small disturbance of specified 
amplitude. The last two aspects of the analysis are mathematically 
modelled and the resulting analytical solutions used to establish 
the strength and distribution of heat sources, for the former, and 
the effect of viscous damping, resonance frequencies, tension, etc., 
on a vibrating ribbon, for the latter. Where feasible, experimental 
checks are made to verify that the designed plate assembly operates 


as required. 
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CHAPTER 1 


INTRODUCTION 


Since the beginning of the twentieth century, much interest 
has been expressed and effort expended in the study of conditions 
affecting the stability of laminar boundary layers as a possible 
explanation for the onset of boundary layer turbulence (see refer- 
ence (1) for a comprehensive bibliography). One such condition = 
the possible stabilizing effect of heating due to a liquid's vis- 
cosity-temperature dependence - provides the motiívation for this 
investigation. Choosing a relatively simple, well-defined experi- 
mental configuration - that of a uniform, parallel, water flow 
over a heated, constant-temperature, sharp-edged, flat plate - 
the significance of this condition is to be examined. It is anti- 
cipated that the empirical results obtained from experiments using 
the apparatus described below might be reasonably compared with the 
analytical solution obtained by Wazzan, Okamura, and Smith (2,3). 

The ensuing discussion, then, outlines the various facets of 
design, construction, and where feasible, testing involved in the 
model simulation. Since a stability study is to be conducted, not 
only should the conditions delineated above for the model be dupli- 
cated as closely as possible, but all those masking the effect of 


the heated flat plate must be reduced to a level of insignificance. 





For the latter, such conditions included the effects of surface 
roughness and waviness of the plate, the effect of a streamwise 
pressure gradient, and the possible importance of buoyancy effects. 
Eventually each of these influenced design considerations, as did 
the more practical requirements concerning the selection of the 
proper materials, ease of access to the parts of the final assembly, 
replaceability of components, elimination of any extraneous vibra- 
tions, maximum utilization of the available testing length, and the 
paramount requirement that the assembly be able to seal against 
water leakage. 

Specifically, the design aspects for the experimental model 
were allocated to three basic areas: 

1. the simulation of a flat plate at zero incidence to a 
uniform, parallel flow with very low turbulence intensity and zero 
pressure gradient; 

2. the construction of a flat plate that would give constant 
temperature on the working surface (the portion of the plate adja- 
cent to the flow over which the boundary layer will build and 
measurements will be taken) with convective heat transfer along its 
length; 

3. the design of a device for introducing a disturbance into 
Ae boundary layer of the plate. Note that no restriction is placed 
on the amplitude of the disturbance, so that both finite oscillat- 
ing waves and the "infinitesimal Tollmein-Schlichting" ones can be 


generated. 





By presenting several possible methods for simulating a 
given condition, and then discussing the criteria governing 
the selection of a particular one, insight is given to the 
interdependency of the numerous design en In cases 
where discrepancies inherently exist between a proposed com- 
ponent design and the assumed experimental configuration, 
modifications to the former are made to optimize overall 
system performance. Since the systematic search for per- 
formance requirements for the composite apparatus followed a 
logical progression once all the required features were 


identified, the ensuing discussion will follow that same order. 





CHAPTER 2 
SIMULATION OF A FLAT PLATE AT ZERO INCIDENCE TO A 


UNIFORM, PARALLEL FLOW WITH LOW TURBULENCE LEVEL 
AND ZERO PRESSURE GRADIENT 


2.1 Testing Facility 
A. Physical characteristics and operating capabilities 

The experimental stability investigation outlined in Chapter 1 
is to be conducted in the 6" x 9", rigid-walled, constant-area 
test section of the closed-loop, low-turbulence water tunnel at 
Case Western Reserve University. Designed and constructed Bee 
Domholdt (4), with a 25-to-l contraction ratio nozzle located 
downstream of a series of damping screens, the facility is capable 
of operating at test section free stream speeds of up to forty 
feet per second with longitudinal turbulence intensities levels 
ranging from .04% at 5 fps to .15% at 35 fps. With a utilizable 
testing length of ade, this reflects a maximum Ke, = 3147. 


4 
The anodized aluminum flat plate (9' wide, El thick, and 


8 
127" long) over which measurements are to be taken is fitted into 
a plexiglass frame mounted on top of the test section so that its 
sharp leading edge extends into the flow approximately en, 
Figure 2.1 shows the location of the plate and the portion of the 
tunnel circuit germane to this discussion. Further explanation 


for the choice of this particular experimental configuration con- 


stitutes the remainder of this paper. 
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B. Uniformity and turbulence intensity of the flow 

Before proceeding with a discussion of the plate itself, it 
would seem feasible to examine briefly the attainability and 
ramifications of the requirement for a uniform, parallel flow with 
a low turbulence intensity. 

Due to the symmetry of the nozzle section immediately preceding 
the test section in the water tunnel, and the fact that the leading 
edge of the flat plate would be located approximately =," downstream 
of the nozzle-test section juncture, it is expected that any second- 
ary "cross'' flows arising from the curvature of the nozzle would be 
insignificant, and the flow would approach the location of the 
leading edge in a uniform, parallel fashion. This expectation is 
confirmed by Domholdt (4) and Nice (5), whose pitot tube and ane- 
mometry equipment measurements, respectively, verified an essentially 
slug-like flow across the test section for all speeds, with the 
expected boundary layer retardation at the sides. The former reports 
a uniformity of to within + 0.5% over the entire speed range, but 
notes that in most cases, better than + 0.25% was observed (see 
figures 6-10 of reference (4)). 

Concerning the free stream turbulence intensity, Schlichting (6) 
points out that based on the results of Schubauer and Skramstad (7), 
unless the turbulence level of the flow is sufficiently small, transi- 
tion will not be caused by the amplification of artificial oscilla- 


tions of certain frequencies, but by the random disturbances already 





present in the flow. Specifically, for a turbulence intensity 
above .03%, transition is caused directly by random disturbances 
rather than by the selective amplification of sinusoidal oscilla- 
tions. If on the other hand, the flow turbulence intensity is low, 
then amplified artificially introduced sinusoidal disturbances of 
a definite frequency can be detected even in the presence of natural 
oscillations inherent in the flow. Schlichting (6), using the 
results of Granville (8), has shown that by combining the experi- 
mental results of Schubauer and Skramstad, Hall and Hislop, and 
Dryden, the difference between the Reynolds numbers at the points 
of transition and instability (plotted as a function of the tur- 
bulence intensity) decreases with increased turbulence, but does 
not equal zero until the turbulence intensity is approximately 2% 
E07 3%. 

Since the establishment of low turbulence levels is so 
important for stability measurements, an attempt to reduce the level 
of the testing facility still further has been made by the addi- 
tion of four more damping screens to the original seven, but 
little or no decrease was observed (5). Apparently, as discussed 
by Schubauer and Skramstad (7), the noise generated as a result 
of the increased power needed to force the flow through the additional 
ns at the same flow rate offset any reduction in free stream 
turbulence intensity. It should be noted that the turbulence inten- 


sity measurements of Domholdt (4) and Nice (5) are only one- 





dimensional ones made in the direction of greatest fluctuation 
(longitudinally) and then extrapolated to three-dimensional, and as 
such, should be considered as upper bounds on the actual free stream 


turbulence intensities. 


2.2 Orientation of the Plate in the Test Section 
A. General location of the plate 

Now that it is determined that a uniform, parallel, low- 
turbulence intensity flow can be established in the test section 
(section 2.1B), the simulation of a flat plate at zero incidence 
to this flow is examined. Logically, the first problem encountered 
concerned the general orientation of the flat plate in the test sec- 
tion, for this ultimately must determine the plate design. 

In an attempt to utilize the existing test section - one in 
which no provision had been built to compensate for the favorable 
pressure gradient that would result from boundary layer build-up 
(e.g., moveable test section walls to change the downstream cross- 
sectional area, boundary layer suction along the walls, etc.) - 
the plate was not put in the center of the test section, contrary 
to the majority of wind tunnel experiments. To do so would cause 
the boundary layer growth on both sides of the plate as well as 
the test section walls, thereby decreasing the already limited effective 
potential flow area and so increasing even more the streamwise 
pressure gradient. As already noted, most of existing experimenta- 


tion for the type problem considered here has been carried out in 





a gaseous medium using a wind tunnel. For a water tunnel where 
sealing and leakage problems arise, as much as possible of the 
associated equipment for the plate normally placed in the tunnel 
on the non-testing side of the plate now must be on the outside. 
Therefore, it is felt that ease of direct access to at least 

one side of the plate would be advantageous. 

A second possibility for the placement was to put the 
plate on the bottom of the test section. However, since it was 
anticipated that the water tunnel should at times be operated 
at low flow rates, and that stability measurements would be 
taken over the heated flat plate, the Rayleigh instability 
problems (where a top-heavy fluid is statically stable under the 
joint influence of conduction and viscosity until the vertical 
temperature gradient becomes sufficiently large) normally associated 
with free convection flow might become significant. Of course, 
similar reasoning followed with respect to locating the flat plate 
on either of the side walls. 

Therefore, the flat plate was installed on the top of the 
test section. In this way, external access was granted to the 
heating elements, pressure taps, thermocouples imbedded in the 
plate, and parts of the perturbing mechanism. 

B. Amount of plate insertion into the test section 
This choice of position for installation, however, created 


other problems when trying to simulate the actual flow phenomenon 





10 
used for the model - problems not present in most wind tunnel 
experiments as already noted. For example, in the latter facility 
where a plate could be placed so that free-stream parallel streamlines 
would SL the leading edge, and test section walls could be 
adjusted to accommodate the pressure gradient, the flow could be 
channelled around the plate so as to simulate a flat plate at 
zero incidence to a uniform, parallel flow. 

On the other hand, for this analysis, so that the boundary 
layer developing in the nozzle section of the water tunnel did 
not interfere with that on the plate, the requirement was im- 
posed that the plate be inserted into the test section to a 
distance greater than the anticipated nozzle-section boundary 
layer thickness. The undetermined quantity, then, was how great 
this distance should be to ensure that the leading edge was at 
least into the region of uniform or "slug-like" flow. According 
to Domhaldt (4), although the boundary layer thickness (based on 
a 2% velocity defect) varied considerably at test section speeds 
less than 15 feet per second (fps), it remained essentially con- 
stant (and greater than for the lower speed range) at about 0.14 
inches for test section speeds above 15 fps (see figure 11, 
reference (4)). His explanation for this observation was that 
the boundary layer remained laminar up to a speed of at least 
7 fps (where $ oC get ); but beyond about 15 fps, the point of 


transition apparently moved farther upstream than it would have 





11 


for a flat plate for the same velocity increase. Consequently, the 
boundary layer thickness remained essentially constant instead of | 
thinning as it would for a fixed position on a flat plate (é £ OE. 
Therefore, in a compromise between assuring that all the pre- 
existing boundary layer could be removed at the leading edge of 

the plate (in a manner discussed in the next section), and the de- 
sire to remove the least amount of water at the leading edge, the 
plate was initially extended into the test section 3/16" or 
approximately .048" greater than the maximum expected nozzle-exit, 
boundary layer thickness. It should be noted that the design of 

the plate assembly allowed for increasing this distance by inser- 
tion of different size gaskets between the plate and plexiglass 
frame holding it in place (see figure 2.2), if later testing deemed 
it necessary. Eventually this distance was increased to approxi- 


mately 9/16" in an attempt to reduce the test section streamwise 


pressure gradient, as discussed in section 2.4C. 


2.3 Leading Edge Suction Requirements 
A. Pumping and piping considerations 


With the irregular constriction now presented to the test sec- 
tion flow by the protruding flat plate, some provision had to be 
made to ensure that the flow remained uniform and parallel as 
before. The obvious solution was to provide just the right amount 


of suction at the leading edge of the plate so that all the flow 
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between the bottom of the plate and the top of the test section 
(including especially the boundary layer build-up from the nozzle) 
would be channelled off, thereby permitting only uniform, parallel 
flow on the testing side of the plate at the leading edge. Hence, 
streamlines in the plane of the plate would split its leading edge 
and a Blasius type boundary layer would form. Note that "suction" 
is used here to denote only the process for ensuring that all the 
flow approaching the plate between the top of the test section and 
the leading edge be channelled into the "suction channel," and 
should not be confused with the term "boundary layer suction" in its 
normal context. By applying suction with the protruding plate, it 
is anticipated that any irregularities in the boundary layer 
(e.g., caused by the discontinuous nozzle-test section juncture 

or a dye injection device to be discussed in the next section) 

can be drawn off also. 

Similar attempts at this type of flat plate design (9,10) 
have resulted in a limitation on the maximum attainable leading 
edge suction, thereby restricting useable test section operating 
speed to less than ten fps, or about one-quarter of the maximum 
capacity of the tunnel. It was felt that this limitation was a 
direct result of inadequate pumping facilities for handling the 
volume of flow required at the net positive suction head (NPSH) 


that existed at the greater speeds. Perhaps this particular 
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aspect of the design required more attention than might otherwise 
have been assumed. 

The whole design and construction of the leading-edge suction 
apparatus was based upon making the head loss through that system 
as small as possible in order to be able to use the smallest (and 
consequently, the least expensive) pump. Since the large volume 
of water that had to be withdrawn made it impractical to discard 
the deionized water of the main tunnel after drawing it off, it 
was decided that this water should be reinserted into the main 
tunnel loop. In this way, the suction apparatus would just form 
an auxiliary loop of the main tunnel flow. It was required that 
the water be reinjected in the tunnel at a position where the 
Static pressure was relatively low and where the injected flow would 
not affect the flow through the test section. Examination of 
figure 13 of reference (4) revealed that the lowest static pressure 
in the tunnel during operation, excluding the test section itself, 
was in the transition section just upstream of the main pump. 

It was anticipated that this auxiliary loop should have 
between the test and transition sections a pump to circulate the 
water (essentially just overcoming the head losses in the system), 
one or more valves to regulate the flow through the piping, and 
some sort of device to measure the volume of water which would 


cover the full range of flows needed. Here again, the selection 
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of the valves and piping size represented a compromise between 

the higher head losses at large flow rates through smaller piping, 
and higher costs and installation inconvenience of larger size 
piping and fittings. After some calculation, it was determined 
that a pipe size of 2-1/2" with associated fittings and the 
physical location of each component would result in an NPSH of 
about 8.7' and require a pump of at least 4 BHP. With these 
figures in mind, it became obvious why difficulty had previously 
been experienced in obtaining sufficient suction for tunnel 
operation at speeds above ten fps (where two 3/4 HP pumps were 
used in parallel). The size of the pump was chosen large enough 
to allow for the possible extension of the plate further into 

the test section (and hence greater flow rate pumping requirement) 
for reasons discussed previously. Consequently, a Goulds 

Model 3643, 2-1/2 x 2-1/2, close coupled, single stage, enclosed 
impeller, centrifugal pump rated at 5 HP at 3500 RPM was chosen. 
It should be noted that a centrifugal pump was chosen rather than 
other types for the following reasons: quiet operation; relatively 
low initial cost and inexpensive and easy maintenance; small size 
for associated output; no excessive pressures even with the 
discharge valve closed: smooth non-pulsating flow: simple construc- 
tion with relatively few moving parts and therefore less chance 

of failure; and high efficiency for relatively large volumes of 


liquid at small heads. 
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Since the size of the valve required for this size piping would 
hardly provide the fine adjustment needed to regulate the amount of 
suction needed at the leading edge of the plate, it was decided that 
a fine control valve would be needed, so a 1-1/2" needle valve was 
placed in parallel with a 3" globe valve. 

The final design of the system appears schematically in 
figure 2.1. All piping with the exception of that on the branch with 
the needle valve is of 2-1/2" polyvinyl chloride (PVC) (the latter 
being 1-1/2" PVC piping). Additional considerations in the place- 
ment of valves included the requirements that the system not have 
to be completely dismantled when the water tunnel has to be partially 
disassembled (to check the screens, remove the test section, etc.), 
couplings be strategically placed so as to allow easy maintenance 
and disassembly when necessary, the minimum amount of piping be used 
to reduce the head loss, and the valves and rotameter be placed 
so as to be easily accessible from one central location - the 
operator's stand. 

B. Attainment of the proper amount of suction 

Due to the wide variation of flows that would be measured for 
the full range of tunnel operation, it was necessary to devise a 
means for measuring approximately how much liquid was being pumped 
sag the suction system. Because the flow in the tunnel could 


be measured but not exactly reproduced from one run to the next, 
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the exact suction needed for any given experiment could not be 
obtained by just reading a dial or gauge in the suction line. 

What was needed was a device that would give reproduceable readings 
for a given suction flow rate which could then be correlated with the 
appropriate, approximate tunnel speeds. Fine adjustments could be 
obtained in other ways as necessary. At this point, three methods 
were examined for determining this flow rate: an orifice meter, a 
venturi meter, and a rotameter or flowrator. The first method 

was discarded as introducing an unwelcomed large pressure drop 
(i.e., a large head loss), and the second was discarded because of 
the limited field that could be covered with sufficient accuracy. 
The third method might have been rejected on the same grounds, 
except that for the rotameter, various size floats could be made 

to cover the flow rates of interest, provided the bore of the 
instrument was sufficiently large. Most importantly, the head loss 
was the least of the three and negligible in all cases. Although 
this last method was perhaps the most costly ordinarily, a Fischer 
& Porter flowrator with a 3'' bore was obtained quite inexpensively 
as government surplus ~ a 1948 NACA item now considered obsolete. 
This instrument was traced by its serial number to obtain the 
specification sheet which indicated that it had been built and 
floats designed for use with gasoline and kerosene. Of course the 
calibration curves, even if they were available, would be of little 


use since it was now intended that the instrument be used in water. 
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As the end result, to cover the upper ranges of the expected flow 
rates through the suction system, it was necessary to design and 
build one float in addition to the one supplied with the rota- 
meter. Calibration for both floats was DELEON by running 
water through the rotameter at a given rate, noting the meter 
reading and the time at which a certain weight of water was con- 
tained in a large cylindrical drum. Note that this weight was 
greater than that accumulated in the time for reaching a steady 
flow. The water was allowed to continue flowing at that same 
rate until a specified greater weight was in the drum, and the 
time between reaching these two weights was recorded. By using 
the elapsed time and the weight of water accumulated, the flow rate 
in gallons per minute was obtained. It should be noted that, as 
expected, these data plot as a straight line for both floats. 

It was felt that although the results obtained represented the 
flow rate to within a few gallons per minute, the curves indicated 
the more important result that the readings were repeatable, 

the primary requirement for this device. As long as a given flow 
rate passing through the rotameter could be measured consistently, 
then for a piven approximate tunnel speed, the flow through the 
suction system could be set near the proper value by observing 


the rotameter reading. 
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With a rough estimate of the proper amount of suction obtained 
by using the rotameter, further refinements were now sought which 
would indicate the exact suction flow rate needed. The importance 
of satisfying this condition as closely as possible was stipulated 
by both matching with the theoretical model to establish the 
Blasius type boundary layer velocity profile, and producing as 
nearly as possible the zero pressure gradient assumed in the 
analysis. To determine the exact suction needed, three methods 
appeared as likely candidates: the injection of dye at the 
leading edge of the plate; a transverse trace made with hot film 
anemometry equipment; and finally, the distribution of pressures 
down the length of the plate. 

The first method, the introduction of dye into the flow (an 
upstream distance of 1.25" between the leading edge and the point 
of injection was used) entailed injection on the same streamline 
that would "split" the leading edge if the flow were unaffected 
by the test section walls, observing visually whether it was 
drawn into the suction channel (too much suction) or passed under 
the plate (too little suction). Since it was assumed that the 
flow entered the test section parallelly, this streamline was 
expected to lie in the same horizontal plane as the leading edge; 
that is, the injection device should discharge the dye at the 


same level as the leading edge. The device for accomplishing this 
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was required to be placed so as not to generate turbulence in the 
free stream. It was believedthat the design of Mateyka (9) did 
not take this factor into consideration. Additionally, it was 
required that this device be retractable from the leading edge 
when the suction system was not in use or when the proper amount 
of suction had been established. Both of these requirements were 
satisfied by the design shown in figure 2.3. Note that by 
inserting the device from the top of the nozzle section, any 
irregularities caused by the flow past the shaft would be drawn 
off through the suction system. 

In the second method, the trace made by the hot film 
anemometry equipment would indicate the boundary layer velocity 
profile. Correlation of this profile with the desired Blasius 
type profile would indicate whether the proper amount of suction 
was being used. It should be noted that this was never expected 
to match exactly due to the inability to obtain a "perfect" zero 
pressure gradient. 

The third method required the installation of pressure taps 
down the length of the plate so that the pressure gradient could 
be monitored. Figure 12 of reference (10) indicated the erratic 
pressure distribution without suction and compared this with the 
distribution attained with the proper amount of suction at the 


leading edge of the plate. Based on these results, it was thought 
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Figure 2.3 Dye injection device 





22 


that by regulating the suction so as to obtain the ee boundary 
layer, the pressure distribution should appear to be approximately 
constant. To monitor these pressures, a series of seven pressure 
taps were installed in the plate, the first being 2-5/8" from the 
leading edge, with each succeeding one 1-1/2" downstream, 

alternately staggered 1-1/2" on either side of the centerline. 

Care was taken in producing the holes to ensure that the edges 

on the working side of the plate were made as squarely as possible 
by drilling from the reverse side and then lapping the plate. At 
first concern developed as to whether the effect of these edges 

not being perfectly Square would introduce any error into the system; 
for certainly, if this error were significant, then the structural 
imperfections from one hole to the next might give erroneous 

readings and explain some of the variation noted by Mateyka (9). 
Franklin and Wallace (11) in evaluating the static hole error for 
holes about the size that were needed here, explained that even 

with a hole with perfectly square edges, a residual error remains 

in the readings due to the local change in the boundary conditions 
and the consequential disturbance of the boundary layer: and so, they 
obtained a plot of (error in pressure reading wall/shear stress) 


| Ur d 
ET) versus the Reynolds number (Re = = , where U. 





T 
is the friction velocity = lic ‚ d is the diameter 


of the hole, and Y is the fluid kinematic viscosity). To obtain 
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an indication of the size of the error for this problem, the shear 


stress based on a Blasius velocity profile was used (6): i.e., 


ie MER: o | Ye 


Therefore, for a free stream velocity of 40 fps at 60°F, a hole 
diameter of .040" and x = 2-5/8", the distance to the first probe 
(these being the conditions that would introduce the greatest error 
for the proposed plate design), it was found that, 


ki = 216,4 and Ty > 8.44 x 107 Ib 





From the (€ / T,, )- Re plot, an error of approximately €= -005904 psi 
orm. 1634 in H,O was found. At this point it should be noted that 

the above curve was based on data obtained with holes for which 
unusual care (using a special technique) had been taken to ensure 
Squareness of the edges. It was further observed by Franklin and 
Wallace (11) in a literature survey on the subject, that hole 
imperfections resulting from the method of hole construction for 

this plate, would increase the error up to as much as five times 

(i.e., up to € = ,0295 psi or .8171 in 1,0). Upon nondimensional- 
izing this latter error in pressure reading with the dynamic pressure, 


a J? 
+e 2 ‚it was found that SEN E = „002740. By com- 


parison, Mateyka (9) plotted the ratio of static to dynamic pressure 
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and reported an uncertainty of + .005. Note that the error found 
here was the maximum expected for any flow configuration. It was 
felt that, since all the holes were identically constructed and 
only the relative pressure readings on the plate were of interest, 
the static hole error could be neglected. 

For the plate designed in this study, a comparison is made 
between the pressure distributions existing with and without 
suction (for, as an example, U: 10.2 fps and A = .701°) is 
presented graphically in figure 2.4, and is discussed more com- 
pletely in section 2.4C. Basically it can be noted that variation 
of the suction rate does not significantly alter the pressure 
distribution, meaning that the pressure measurements will not 
provide an accurate indication of the exact suction needed for 
a given free stream velocity. 

In fact, evaluating the methods for "fine tuning" the suction 
rate discussed above, the following order of utilization evolves - 
from the least to the most accurate: 

l. pressure measurements 
2. dye injection 


3. boundary layer traverses with 
anemometry equipment 
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Note: Successive taps alternate 15" either side 
of centerline. 
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Figure 2.4 Effect of leading edge suction on the flat plate 
pressure distribution for A= .701° and U = 10.2 fps. 
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2.4 Pressure Gradient 
A. Effect on flow stability 

With its fixed walls, the one aspect that the constant-cross- 
sectional-area test section cannot compensate for is the favorable 
streamwise pressure gradient occurring naturally because of 
boundary layer growth, and augmented by plate insertion into the flow 
(see section 2.2A). Therefore, since the proposed experimental 
configuration assumes a zero pressure gradient, before implementing 
modifications to the existing facility, it is essential that the 
effect of a pressure gradient be examined. 

The strong dependence on the existing pressure gradient of 
the limits of stability and the degree to which a small disturbance 
would be amplified have been shown experimentally by Schubauer and 
Skramstad (7), revealing that decreasing downstream pressures cause 
a damping of velocity fluctuations in the laminar boundary layer. 
Earlier, Schlichting (12) theoretically analyzed for stability, 
the velocity profiles in convergent and divergent channels resulting 
from the potential flow given by U(x) = ws - (constant) X 
with a similar conclusion. Later, Schlichting and Ulrich (13) 
carried out stability calculations using a modified Pohlhausen 
approximate method and determined the critical Reynolds number as 
a function of a dimensionless shape factor, defined as: 


sl) 
De 


Y 
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(where for the purposes of the ensuing discussion, 5 is the 
maximum boundary layer thickness on the plate for a given luis: ; 
occurring, of course, at the trailing edge of the plate). This 
parameter can be interpreted physically as a ratio of the pressure 
to viscous forces. 


It can be seen from a plot of Schlichting and Ulrich (13) 





x 
of the critical Reynolds number Res w = Ze versus the 
2 cr} 
shape factor Ay = 3 au , the small change of A away 
X 


from the zero pressure gradient case for this favorable one 
results in a correspondingly small increase in the critical 
Reynolds number. For example, if N = „23, an apparently 
characteristic value by subsequent theoretical computations (see 
section 2.4B), the increase in the critical Reynolds number is 
approximately 90. From the solution of Wazzan, Okamura, and 
Smith (3), such an increase could be expected from only a 
AT = er - To = 3°F (see figure 3.1). Since it is expected 
that experimental observations will be made at temperature dif- 
ferences greater than this (see section 3.2A), the effects of 
heating the plate should cause a more significant increase of the 
critical Reynolds number than the existing favorable pressure 
gradient. 
B. Design estimate of the magnitude of the pressure gradient 

For the particular case involved here, with fixed walls, 


boundary layer growth on the sides of the test section would 
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accelerate the free stream flow as if it were a perfect fluid 
in a convergent channel. Since the thickness of the layer (and 
consequential decrease of the effective free stream area) is 
dependent on the free stream velocity, one would expect 
that the tunnel speed (in addition to other factors such as 
the temperature of the flow, geometric configuration of the 
test section - it does naturally converge slightly with a 
difference between the inlet and outlet areas of A, = A = 
0.067 A etc.) would affect the magnitude of the pressure 
gradient. 

Making the following assumptions, an indication of the 
effect of the pressure gradient, the shape factor IN , can be 
calculated approximately: 


1. A Blasius laminar boundary layer build-up on the surface 


of the plate which implies (6) 


ga = 


& (x) = 5.0] YE 


o0 


x LE 

Se 72 5 

we 

2. A fully developed turbulent boundary layer on the three 
remaining test section walls (which are treated as flat plates, 


i.e., corner effects occurring due to finite size of test section 
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are neglected) where the boundary layer thickness Š , based on 
a 2Z velocity defect, increases from approximately 0.14 inches 
at the test section inlet (for Ua = 15 fps) to .25 inches at 
the outlet (4). It is further assumed that the values measured 
at the nozzle exit (i.e., at the test section inlet) are 
Pee roximately the same as those at the leading edge of the plate 
(the point on the plate which encounters the free stream first, 
and from which all measurements relating to the boundary layer 
are to be measured in a streamwise direction), thereby giving 
a conservative value for A . 

3. Turbulent boundary layers are defined by the 1/7-th 


power velocity distribution law (25) where 


from which can be determined 5. = 8/8 ° 

4. The potential velocity gradient is linear. 

5. The plate is parallel to the top of the test section 
and inserted to a distance of 3/16". 

Since the displacement thickness 5” , Ís defined as the 
amount by which the external potential flow is displaced by the 
decrease of the velocity in the boundary layer, and since the 


potential flow distribution essentially determines the test 
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section pressure gradient, by calculating the difference in 
"effective cross-sectional areas" (caused by the boundary layer 
build-up and the natural convergence of the test section), the 
potential velocity gradient can be determined. Calculation of the 
shape factor follows directly. 
For example, for the unheated flat plate at a temperature 
of To, = 60°F, 
Veo (458) ¿a 
+0 Ze 
29 We, 
Since there is a laminar to turbulent boundary layer transi- 
tion at the nozzle exit between the tunnel speeds of 10 fps and 
15 fps (4), it is difficult to ascertain exactly what the boundary 
layer thickness on the walls of the test section will be at the 
leading edge of the plate, due to the finite distance between that 
point and the nozzle exit (where the boundary layer thickness is 
known). Based on previous assumptions this should be no greater 
than 0.14 inches. Assuming that this value may be used for a 
tunnel speed of 10 fps, the shape factor may conservatively be 


calculated. 


U (fps) ae 
LO EL 
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Notice that the values for A are quite consistent, under the 
assumptions that have been made here, with the experimentally 
determined maximum value of Mateyka (9), EN = „233. Note 
however that this value is based on his measurements, conducted 
in the same test section, where the highest free stream velocity 
used was oe = 7.63 fps. 


C. Experimental measurements and procedures for reduction of 


the pressure gradient 


Upon initial installation of the plate assembly into the tunnel 
(i.e., at zero angle of incidence), it was found that the actual 
favorable pressure gradient in the test section was apparently 
significantly larger than that previously estimated. Possible 
explanations were advanced: (1) insufficient suction, but it is 
found that once suction is applied and the stagnation region 
removed from the leading edge, changes in the suction rate do not 
significantly affect the pressure gradient (see figure 2.4 for 
example); (2) misalignment of the working surface of the plate 
relative to the bottom of the test section, but measurements later 
indicated only slight convergence, and that no worse than the area 
change already noted in the preceding section; (3) boundary 
layer build-up in corners which was not adequately accounted for 
in the analytical estimation of section 2.2B; (4) consideration 
of the pressure reading obtained from the last downstream pressure 
tap, yielding misleading results for the pressure distribution over 


the major portion of the plate. 
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The significant deviation of the last downstream pressure tap 
reading (as seen in figures 2.4, 2.5, and 2.6) is attributed to the 
close proximity of that tap to the trailing edge of the plate 
(.375"). It is felt that the curvature of the stream as it seeks 
to flow around the square edge of the plate'’s downstream end causes 
a lowered static pressure in that region. Hence, although this tap 
does not provide data consistent with the essentially linear pressure 
gradients obtained from the remaining locations, the data so 
obtained does show that subsequent adjustment of the plate's angle 
of incidence fails to correct the large pressure drop, and that 
small pressure gradient data must be taken near the leading edge 
or central portion of the plate so as to remain away from the large 
drop. Therefore, comparing experimentally measured with the 
analytically/empirically calculated values of the shape factor, 

AN , presented in section 2.4B, the former is found to be 
larger but essentially constant for all free stream velocities. 
Without considering the last point, on an average, aS = .8977 
whereas consideration of the last tap reading gives an average 
value of IN = 1.3616, an increase of approximately 522. 

Of the four explanations advanced, the last two seem to be 
the most likely for the disparity between the calculated and 
Be rved pressure gradient. 

At this point, several possibilities were considered for 


reducing the test section streamwise pressure gradient: 
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Figure 2.5 Flat plate pressure distribution for zero angle of 
incidence to the mean flow at various flow velo- 
cities, 
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Figure 2.6 Flat plate pressure distribution for different 
angles of incidence to the mean flow at U = 8.0 fps. 
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(1) modification of the test section. Fix within the test 
section a "flexible false bottom" which would extend from within the 
nozzle section, through the test section to the split diffuser 
section, the contour of which would be adjusted by deformation at 
several positions along the length. In this manner, the adjustable 
floor contour would compensate for the decreased effective area for 
each tunnel speed. 

(2) inclination of the plate. Utilizing the designed versatility 
of the plate assembly, to incline the working surface to the mean 
flow without drastically altering the external piping, sealing of the 
plexiglass frame supporting the test plate, etc., the original 
1/32" cork gasket sealing between the plate and the plexiglass 
frame would be replaced by a thicker, pliant rubber gasket material 
(3/8" medium density closed-cell neoprene sponge rubber eventually 
was chosen) (see figures 2.2 and 2.7). In this way, the angle of 
the plate could be adjusted by merely compressing the gasket more at 
one end than at the other while the frame and piping remained fixed 
in place. 

Either method of correction would be made to compensate for 
boundary layer growth on all four sides of the test section - not 
just that on the plate - and so the required correction would be 
much larger than that of only the displacement thickness of the 


boundary layer on the plate. 





36 


Although it was not anticipated that the second correction proce- 
dure would completely eliminate the pressure gradient as perhaps the 
first might have, it was thought that the reduction would be signi- 
ficant enough to give it preference over the more mechanically involved 
modification of the test section. It was not expected that the 
small angle of incidence would significantly affect the results. 

An additional advantage claimed by the second procedure can be 
utilized; it would require that the plate be inserted further into 
the free stream than the 3/16" calculated in section 2.2B. Measure- 
ments of SES intensities indicate that the effect of 
this turbulence is discernable further into the free stream than 
the measured boundary layer thickness (5); the additional distance 
between the leading edge and the top of the test section gained 
by the thicker gasket would ensure that the "turbulent intensity 
boundary layer" could be drawn off as well as the velocity boundary 
layer. This point will be clarified in subsequent discussion. 

Since the experimental analysis herein is directed toward 
confirming that the designed apparatus functions to duplicate a 
specified model, attention is focused primarily upon establishing 
trends - that is, the qualitative rather than quantitative behavior 
of the equipment. Therefore, no excessive effort was expended 
to amass a huge collection of data for different suction rates, 


angles of plate inclination to the flow, etc. Instead, the plate 
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assembly was operated in various configurations (see for 
example, figures 2.4, 2.5, and 2.6) to see how assis it could 
satisfy the required conditions, then assessing the significance 
of any disparity between the two. To this end, the following 
observations are made. 

There are two factors which apparently cause the appearance 
of an unfavorable pressure gradient in the test section: (1) no 
suction at the plate's leading edge, and (2) increased inclina- 
tion of the plate to the flow. As may be verified by figures 2.4 
and 2.6, this gradient is generally small with respect to the 
total head, with the maximum pressure occurring at or slightly 
upstream of the plate's midpoint. 

l. It appears that once suction is provided at the leading 
edge (so that the flow no longer encounters the stagnation region 
there as if for a no-suction configuration), wide variations of the 
suction rate affect the pressure gradient only slightly (see 
figure 2.4). 

2. As displayed in figure 2.6, although an increased plate 
angle of incidence to the flow does tend to decrease the pressure 
gradient initially, further inclination after a certain angle will 
introduce an adverse pressure gradient and apparently not decrease 


the magnitude of the overall gradient. 
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At increased free stream velocities, the effect of Us 
on the pressure gradient seems to decrease; that is, analysis of 
figure 2.5 shows that the value of Te JZ decreases at a 
given location on the plate as o increases. 

Also, although trends in the pressure distribution were 
apparent, appearing consistently in the limited amount of data 
taken, within the limits of the accuracy of the readings, the 


pressure gradients for a free stream velocity less than five 


feet per second can be taken to be zero. 


2.5 Design Considerations for the Plate 
A. General discussion 

The details of design for the plate were purposely kept as 
simple as possible to facilitate construction and to allow versatility 
(see figure 2.7 for ensuing discussion). Guided by the desire that 
the components of the plate assembly be replaceable without having 
to rebuild the whole assembly, the final design consisted of 
three parts: the actual flat plate, a plexiglass supporting frame, 
and a transition section which would channel the flow from a 
rectangular cross-sectional suction channel to a standard size 
pipe fitting. 

The plate would fit into a slot in the frame, would be held in 
place by a series of studs embedded in the plate, and would be sealed 


with a flat gasket "sandwiched" between the two. Pressure tap 
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Figure 2.7 Exploded view of the parts of the plate assembly 
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fittings also embedded in the plate would pass through holes in the 
frame. All such holes were made oversized and studs, rather 
than bolts, were used in securing the plate to the frame to permit 
the variable longitudinal inclination of the plate relative to 
the frame required for the pressure gradient elimination discussed 
in section 2.4C. In this way, plate pieces could be easily 
changed without disrupting the rest of the apparatus, sealing 
problems would be minimized, access to the pressure taps would be 
unhampered by a top plate (i.e., there would be no problem with 
having to seal around where the pressure taps come out of the plate 
due to the gasket seal). Similarly, for the magnet chamber, access 
to the magnets would be unhampered by a sealing enclosure, and so 
virtually any reasonable size magnet required for the proposed 
electromagnetic perturbing mechanism could be accommodated. The 
dilemma then became a compromise between how far forward the magnet 
chamber could be placed to the leading edge while still allowing 
the suction channel at the leading edge, and the requirement that 
there be sufficient structural strength between the two so that 
failure would not occur when operating at maximum tunnel conditions. 
The whole apparatus had to be constructed so that extraneous 
vibrations of the plate resulting from the large flow rates were 
kept to a minimum. To this end, for example, rather than making the 


leading edge of the plate completely planed off at the angle of 
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ten degrees, small "lips" were left on both sides. Rather than the 
plate fitting flatly on the frame, the latter was grooved so that 
the plate would be held in place not only by the walls of the test 
section, but by the recess in the frame as well. 

In designing the suction channel, a nwmber of factors were 
considered. To utilize the testing length to the fullest, the 
suction channel had to be as far forward as possible, thereby 
dictating ar the flow be drawn off vertically (as opposed to the 
more favorable small angle with the incoming flow) to accommodate 
the perturbing mechanism that would have to be placed on the plate. 
The width of the suction channel was maintained essentially 
constant throughout the turn to ensure that the flow across the 
inlet remained as uniform as possible. To reduce the friction 
losses with the relatively large flow rates that would be experienced 
at the maximum capacity of the water tunnel, the cross-sectional 
area of the suction channel was increased from the inlet (A = 1.69 ioe 
through the plate and transition section, to the piping (A = 4.26 ay 
that carried the flow to the pump. As the flow in the channel would 
initially be turned through a right angle, even a cursory examination 
of the friction losses around a mitred corner indicated the advantage 


of curving the channel in the flow direction. 
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B. Selection of plate material 


With the theoretical design considerations taken care of, there 
remainedthe problems of a more practical nature - the selection of 
the material for the plate and then the Petes of the plate itself. 

For purposes of this investigation, a hard anodized 6061 T6 
wrought aluminum plate was found to satisfy all the requirements for 
the test piece (14 - 20). Specifically, since no large localized 
stresses were anticipated, a softer, more easily machineable 
core material could be used under a hard surface coating. The 
latter, AL, oF , while strongly adhesive to the aluminum 
substrate, provided a highly dielectric, corrosion and abrasion 
resistant surface. Both materials were nonmagnetic to avoid 
shorting the magnetic flux path needed for the disturbance 
mechanism as discussed in chapter 4. Further motivation for 
this type of coating stemmed from its successful application on 
high speed airfoils where the leading edge would be susceptible 
to damage due to the surface attrition of dust and rain. 

The only feature of the anodic coating which affected the 
physical specifications for the plate was the increase in each 
surface dimension by about one-half the thickness of the anodic 
coating produced. Consideration of all the properties of the 
coating was felt essential to ensure that the plate would perform 
properly, and the results of that study are illuminated in 


appendix A. 
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C. Surface roughness 


Since this plate was designed for use in a stability study, 
the effects of surface roughness and waviness on transition had to 
be considered - specifically, the fact that roughness elements 
would generate disturbances in the laminar boundary layer in 
addition to those already present there due to the flow turbulence 
level. If the former disturbances were larger than the latter, 
then a smaller degree of amplification would be required then for 
the smooth plate surface to effect the transition. Conversely, if 
the surface were very smooth, disturbances generated by surface 
roughness elements would be less than those already present due 
to the free stream turbulence, and the presence of roughness 
would have no effect on transition. Therefore, it was necessary 
to determine the degree of surface roughness permissible for the 
flat plate. 

In the limited number of experiments on the effect of distri- 
buted surface roughness on transition, the results have been some- 
what difficult to analyze due to the necessary experimental treat- 
ment of the leading edge. According to Lin (21) in analyzing 
Holstien's study (22), a Reynolds number based on the free stream 
velocity, the nominal roughness height, and the existing viscosity, 
must be less than about 100 + 50 for the surface to be considered 
aerodynamically smooth. The experiments conducted later by Feindt 


(23) indicated that the grain size, K, , in the presence of 


S 
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sand roughness had a significant effect on the critical Reynolds 
number only after the value 


a= eo Ky = 120 


a m 


had been exceeded. Prior to this point there was no apparent 
change in the critical Reynolds number, but after, the value de- 
creased steeply. For this problem, using maximum operating con- 
ditions, Ke > 334 microinches. Hence if the grain size were 
kept less than this value, the distributed surface roughness should 
have no effect on transition. Tests on the surface of the plate 
in its final form (including the surface anodization to be dis- 
cussed later) were conducted, and the nominal average surface 
roughness was found to be 30 microinches, certainly well within 
the limits discussed above. (This would correspond to U. Ks =8,9 
y 


as opposed to the limit of 120 noted above.) 





CHAPTER 3 
THE GENERATION OF HEAT REQUIRED TO MAINTAIN A CONSTANT 
TEMPERATURE FLAT PLATE IN A UNIFORM PARALLEL STREAM 
3.1 Selection of Heating Method 
To model a heated, constant temperature flat plate at zero 
incidence to a parallel laminar flow, a special method for heating 
this plate had to be devised, dependent upon assumptions concern- 
ing the temperature dependence of flow parameters and the spatial 
dependence of thermal properties of the plate (i.e., homogeneity 
of material), the rate of surface heat flux as determined by the 
formation of the thermal boundary layer, the transient effects 
governing the time required to attain steady state, the temperature 
difference between plate surface and free stream and associated 
buoyancy effects, and so forth. Experience and theory have shown 
that, due to the large temperature gradients near the leading 
edge of the plate caused by the thinness of the thermal boundary 
layer, the rate of heat transfer for a given temperature would be 
greatest there, decreasing with distance down the plate as the 
boundary layer became thicker. Therefore, regardless of any 
assumptions made in trying to simulate the phenomenon, if a con- 
stant temperature were to be maintained on the surface of the model, 
spatially varying amounts of heat had to be supplied to the plate 
for transfer to the nel, with the greatest amount being at the 


leading edge. A survey of the literature indicated several possible 
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methods for accomplishing this heating: passing a vaporized liquid 
(gas) or a mixture of the vapor (gas) and liquid on the reverse 
side of the object to be heated or cooled (to utilize the constant 
heat of condensation at the appropriate pressures) either over 

the whole surface or through a series of tubes (much like a heat 
exchanger) (9, 10, 24); the use of a number of rheostat controlled 
strip electric heating elements (25, 26). With each method, there 
were problems in establishing and maintaining the constant tempera- 
ture. For example, for the installation under consideration where 
the plate would be mounted horizontally with the working surface on 
the bottom, heating the plate by the condensation of steam would 
cause the formation of film or drop-wise condensation on the upper 
surface of the plate, thereby affecting the uniformity of the heat 
transfer. When using either the passage of a fluid through a series 
of coils or the electric resistance heating of a series of heating 
elements, the uniform heat flux on the working surface of the plate 
is being approximated by a finite distribution of sources, so the 
temperature can never be constant; that is, there will always be 
some "waviness" in the spatial temperature distribution on the 
plate. However, owing to the success of several experiments (25, 
26) where heating of a flat plate was accomplished by use of 

strip heating elements, and the relative simplicity of the equip- 


ment needed, this method was chosen for this experiment. 





3.2 Implimentation of Chosen Method - Determination of 
Distribution of Heating Elements 


A. Physical limitations affecting design 

With the method of heating now chosen, it remained to actually 
design the plate so that the internal heat sources would closely 
approximate the actual heat transfer. In this way, the spatially 
periodic variation of temperature on the working side of the plate 
caused by the discrete distribution of sources separated by 
adiabatic regions would be made as small as possible within the 
limits of the physical size and capacity of the heaters. For 
this design, a temperature difference, AT = [Tisas T | = 1°F 
along the length of the plate was sought, an improvement on the 
"constant temperature” heated flat plates used previously in this 
water tunnel by Mateyka (9) and Morris (10), who reported a varia- 
tion of approximately 8° - 15°C. Note that in formulating this 
problem, it was not stipulated that both the surface temperature 
and heat flux be specified (the two conditions that should match 
with the exact laminar thermal boundary layer solution) to determine 
the distribution of sources. Since the equations defining conduc- 
tion in the plate were elliptic, it could be shown that the imposi- 
tion of Cauchy conditions on the surface of a closed domain over- 
specified the problem, and for an open domain, made the solution 


unstable (27). 
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Since the shape and size of the plate were predetermined by 
the test section and the requirement that the heated plates be 
interchangeable (i.e., all be able to fit into the frame in the 
test section), the problem was reduced to ide the place- 
ment and strengths of the heaters. Admittedly, the problem might 
be solved with very little computation by building a plate with 
an arbitrary placement of the heaters, and varying the heat out- 
put of each until the surface temperature was as constant as 
possible for that distribution; but it was felt that an analytical 
solution prior to the construction of the plate would permit the 
selection of spacing and distribution of strengths of the sources 
to better simulate a constant temperature, and the "fine tuning" 
could still be performed to remove any minor disparities. 

To this point, it has been assumed that the amount of heat 
required to maintain a free-stream-to-wall temperature difference 
could be supplied, but no indication has been made regarding 
its magnitude. Specifically, the maximum temperature difference 
that must be maintained is required. Since the tunnel is limited 
in velocity and the plate in available testing length (the mini- 
mum being Te position of the perturbing mechanism as measured from 


the leading edge of the plate) to, 
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Figure 3.1 Plot of Re versus wall temperature based on the 
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For an initial estimate of the temperatures required, by 
plotting Re „a versus wall temperature from the data of Wazzan, 
Okamura, and Smith (2) (figure 3.1), and then superimposing these 
limits of Rex , it can be seen that the wall temperature must 
remain between that of the unheated plate and that with a difference 
I IK = 22,5°F. For design purposes, then, a maximum 
temperature difference of 25°F is used. 


B. Analytical solution to the steady-state conduction equation 


 — im ee  —n 


in a rectangular domain with rectanpular heat sources 

Prior to a mathematical analysis of the problem, it was thought 
necessary to have some idea of the design in mind. The one 
selected for this case, the anodized aluminum plate, is the same 
as that shown in figure 2.7 with a series of heating elements 
embedded in it. It was felt that the degree of sophistication in 
mathematically modelling the plate to determine the distribution of 
internal heat sources required for a given heat transfer rate from 
the bottom was optional toacertain degree, but eccentric geo- 
metries, assumptions of inhomogeneity in the plate in either the 
vertical or horizontal direction, and the like each only made the 
solution more complex, but not necessarily more meaningful for 
design purposes. For purposes of this analysis, the following 


assumptions were made: 
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l. The plate is of rectangular shape with length AM, 
thickness h , and width w , which can be treated two-dimension- 
ally longitudinally ( R> and w> h A 

2. The heat would be supplied by a finite series of heating 
elements modelled here by a series of rectangular sources. 

3. The top, front, and back of the plate are assumed to be 
insulated, permitting no flow of heat, and the bottom with the 
anodized coating acts as a partial insulator - a resistance to 
heat flow. The efflux of heat from the bottom can be approximated 
by the Blasius-type constant wall temperature similarity solution 


for constant property boundary layer flow (6); 


3 
3 = 332 k VR (be (To-To) L (3.2.1) 
V | x 

Admittedly. however, the assumption of constant wall-to-free stream 
temperature difference and monotonically decreasing X dependence 
implicit in this formulation would not be compatible with the tem- 
perature difference, O , Obtained. That is, due to the discrete 
distribution of sources, it is expected that some oscillatory 
component must be represented in the solution, the amplitude of 
which would decrease directly as the spacing between these sources. 


Perhaps the heat flux would, then, be better approximated by 





53 


3 
jë = OO aie ÍP- | Ue O 


\ 


1 
s. 
DA 


A T > - 
re 1 q cos Tr 2 C un) 


Kan 


where. 
f(x... o X: > Jin L se) 
represents an unknown amplitude of oscillation and is a function 


of the spacing between successive sources, 


and 
E m 

are step functions restricting a particular sinusoidal contribution 
to a particular interval. Of course, the exact solution would require 
the extensively more involved solution of the coupled conduction and 
convection problem. The goal of this analysis is to reduce the amli- 
tude of these oscillations below a given value, and, ideally, the 
"corrective" term should help in this reduction. For purposes of 
this study, however, the amplitude is assumed to be small (i.e., 

Eee 1), so only the leading term is retained. Modification in 
strength and position of the sources is made until this assumption 


is justifiable. 
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4. The plate would be a composite consisting of aluminum 
of thickness b and thermal conductivity P and an aluminum 
oxide coating of thickness C and thermal conductivity K. ; 
where oC and Ka >> sa : 

5. Each region is homogeneous (i.e., uniform thermal proper- 
ties) but discontinuous at the interface between the two. 

6. There is perfect thermal contact at the interface. 

This model is shown in figure B.l. 

Before proceeding with the solution to the problem or intro- 
ducing the results of sample calculations for various source con- 
figurations, a few observations should be made regarding the 
validity of the expected solution. As noted in assumption 3 of 
this section, the analysis assumes a heat flux varying inversely 
as the square root of X (equation (3.2.1)). Implicit within this 
assumption are (1) a constant wall-to-free-stream temperature 
difference NS , and (2) the requirement that there be 
an infinite amount of heat flux as X AE As a practical 
matter, the latter statement is obviously not realistic. In fact, 
as the leading edge is approached, due to the inability to physi- 
cally provide sufficient source output in that confined space, 
the wall temperature would tend toward the free stream temperature, 
A F: . That is to say, with virtually no temperature 


difference between the two, the heat flux would tend to zero in 
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that region. It is not unexpected that this mathematical and 

physical paradox should appear when keeping in mind that the heat 

flux used for this model is derived from the energy ee layer 

equations which are invalid near the leading edge ( X 7 0). ` 
Stating this conflict in another way, whereas the model 


TT 


requires lim A o?’ physical intuition requires that 
E AR: 

i, ear be finite, thus violating the assumption that the wall- 
KO, yx! 
to-free-stream temperature difference be constant. In short, 
the mathematical solution in the region near X = 0, as shown 
in figures 3.2 and 3.3 for example, should not be taken as being 
physically meaningful. 

Additionally, the "hot" regions existing in the neighborhood 
of the sources and the relatively "cold" regions between them, by 
having a larger or smaller Tur lo , respectively, would índuce 
a corresponding increase or decrease in the heat transfer in those 
regions. Unfortunately, this feature cannot be incorporated into 
the present analysis without the simultaneous solution of the coupled 
problem, as pointed out in assumption 3 of this section. However, 
since it will tend to decrease the amplitude of the oscillations 
between sources, the solution obtained herein must represent the 


"worst case” (i.e., an upper bound on the amplitude of the oscilla- 


tion for a given source distribution). 
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Keeping these limitations in mind, a concise restatement of 
the mathematical model is made, and then the temperature distri- 
bution is sought. Heat is generated by a finite series of 
rectangular heat sources in one region of a two-region, composite, 
rectangular slab which is insulated on three sides, but which 
transfers heat by convection to a laminar boundary layer on the 
remaining side. 

The analytical solution is obtained by first determining 
the proper source function (or since the temperature gradients 
on all sides are the prescribed boundary conditions, the Neumann 
function) which, when found, could be used as a kernel in an 
integral equation incorporating internal sources and boundary 
conditions to define the temperature distribution within the 
Mate (see appendix B, sections B.l and B.2 for complete deriva- 


tion). Denoting the aluminum as region " @ " and the anodized 


coating as region " C ," the Neumann functions, N , are found 
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where 1) Vor š = `] m #0 (3:234) 


2 > Imn 
2) Sources are assumed to be in region ]. 
3) SES en are determined by: 
Cot en P + Ka cot enc = © (222.5) 


Ke 


Assuming that the heating elements could be represented by a series 


of step functions 


Sa (x,y) = 2 S Hl] - 
Hle- STB HL yy] - ALY 
(ye SD] ae 
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where O= strengths of the L° source 


(X; +40) 
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To determine the eigenvalues for the problem, characteristic 


th 


position of the center of the L” source 


width of each element 


thickness of each element 


values were used in the eigenvalue equation Stated above. With 


b= .25 C= 002" 


= Btu 
~ 17 — ao. K = = .317 Xi r os 
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the following eigenvalues were found for En ' Waal cee 


t: 


dar 


1 


° SO 


SAO O rr 


r 


After closer scrutiny, it was realized that En was approximately 

Br A na cena ey the 
real eigenvalues were less than this, but it was felt that the 
increasing difference between the real and approximate values was 
insignificant, due to the relatively rapid convergence of the 
solution with increasing On ° 


Therefore, making the approximatiíon that en” — , the 


source functions became, 
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This significant reduction of the equations implied physically 

that heat transfer through the anodized coating would be essentially 
one-dimensional (normal to the bottom surface), and little or no 
enhancement of longitudinal heat transfer could be expected due to 
the presence of the anodized coating. In short, the plate could 

be treated as a homogeneous material. 


In this new formulation the temperature distribution became 
GO 


2m. A Cos e Cos ei sin mE 
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m=1 
where 3 ( X r. ) is defined by equation (3.2.1). 


C. Asymptotic seríes solution to conduction problem as an 


indication of the required source distribution 


In an attempt to obtain a rough indication of the size and 


position of the sources required, one can examine the heat flow in 
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the plate on a still simpler basis - the assumption of a one-dimen- 
sional heat flow. The justification for such a formulation fol- 


lows. 


The governing equation and boundary conditions for the 


rectangular domain are 
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and S (y=b) is the quantity to be determined. 4 is again 


defined by equation (3.2.1). 


Nondimensionalizing with 
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the system of equations (3.2.10) becomes 
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Further assuming that A can be expanded in the parameter ( 5/4 ), 


A 


aan 


and the following set of ordered equations evolve. 


> O. - O 
1 
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and so on. 


Solving for each order and applying boundary condtions away from 


the end points, it is found that to second order, 
bak S I a 
= + = 
E) de en ] | TM sone 
= N nie, Zn. u 
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Notice that this solution is valid as long as no attempt is 


C2500) 


made to satisfy boundary conditions at either end. This suggests 
the presence of a boundary-layer type solution or a singular per- 
turbation problem, an observation substantiated by Van Dyke's rule 
for the uniformity of a solution in a given domain (28); "A per- 


turbation solution is uniformly valid in the space and time 
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coordinates unless the perturbation quantity is the ratio of two 


' 


lengths." Note further that in solving the problem thus far, 
the basic equations have been changed in nature from elliptic to 
parabolic, thereby permitting the imposition of Cauchy boundary 
conditions on one boundary (see section 3.2A) to see what the 


conditions should be on the opposite one. 


From equation (3.2.1) 
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Evaluating the solution at the top of the plate, N = ], and 
realizing that except very close to the leading edge, Too 
a good representation of the required heat flux needed at the 
top, Q ( X, b ), is glven simply by, 

AA 4 0) 
since (3) << | 
A 

To represent the continuous heat flux distribution by a discrete 
distribution of heating elements, ()( x, 6) is integrated over 
some range, AX , which is determined by the size and output capa- 


city of a heating element. That is 
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Qlos, b)= ctt K, (PSE (Tu-A NE) 


(352210) 


D. Discussion of plate temperatures for different source 
distributions 





To assess the capability of various distributions of heat 
sources for providing a required temperature distribution on the 
plate, two examples are presented as sort of "numerical 
experiments." For both, it is assumed that twelve heaters supply 
the required heat output, and further, that each can individually 
provide any specified amount. Additionally, although three variables 
would ordinarily be available for adjustment (i.e., horizontal 
and vertical spacing, and source strength) due to the limitation on 
the vertical spacing by the thinness of the plate, only the heater 
Strength and horizontal spacing are assumed adjustable. Physical 
parameters applicable to both examples are presented in table 3.1. 

The first example specifies a distribution of constant 
strength heaters with variable spacing (such as used in reference 
(26)). This spacing is determined by finding the total heat flux 
required per heater, and then using the one-dimensional-analysis, 
equation (3.2.14), to estimate the distance along the plate re- 
quired to transfer that much heat by convection. The heater is 


then assumed to be placed in the middle of this interval. Strengths 
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120252 

25, 

1002" 

.340 Btu/hr ft °F 
117 Btu/hr ft °F 

.597 Btu/hr ft °F 
8.03 

20.0 fps 

25,027 

.001 


5 


A sec 


25, 


015" 


95,197.00818 a 
in 


50.9983 e 
in 


Table 3.1 
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and positions of the heat sources for this example are listed in 
table 3.2. Note that the first three heaters are not spaced in 
accordance with equation (3.2.14) due to their width. That is, 
the theory would require that they be spaced successively at 
X = ,0425", .2125", and .5525". With the heater width speci- 

fied as A = .25", such spacing is impossible, so they are 
positioned as closely as physical dimensions will permit. The 
resulting temperature distribution is plotted as a percentage 
of the wall to free stream temperature difference in figure 3.2. 

The second example specifies a distribution of variable 
strength heaters with constant spacing (such as used in reference 
(25)). The required strengths are found in much the same manner 
as for example 1, except that here, a distance along the plate 
is specified and the heat transferred by convection in this 
interval evaluated. The appropriate source, in turn, is assumed 
to supply this quantity of heat and to be positioned in the center 
of the interval. Table 3.3 lists the source strengths and posi- 
tions for this example. Temperature variations over the entire 
length of the plate as well as details of this variation in a 
characteristic interval are seen in figures 3.3 and 3.4, 
respectively. 

Discussing these two examples qualitatively, with reference 


to figures 3.2 and 3.3, it is noted that both formulations reveal 





Source 


10 


11 


2 


66 


watts 
en) 
in 


Strength ( 
7,933.0840 
7,933.0840 
7,933.0840 
7,733.0840 
7, 93320846 
7,933.0840 
7,933.0840 
1,9 3320840 
7,933.0840 
7,933.0840 
7,933.0840 


7,933.0840 


Table 3.2 


x (in) 





1250 


. 3750 


.6250 


1.0630 


1.7440 


2.5950 


3.6150 


4.8060 


6.1680 


7.6990 


9.4000 


11.2720 





Source 





10 


11 


12 


67 


watts 


Strength Co (C 5? 


in 
27,199.0980 
11,266.2350 
8,644. 8866 
7,287.9764 
6,420. 8362 
5 804.8797 
5, 338.1377 
4,968.6176 
4,666.6274 
4,413.8061 
4,198.1029 


4,987.6405 


Table 3.3 
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Figure 3.2 Plate surface temperature distribution for constant 
strength sources with variable spacing 
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Figure 3.3 Plate surface temperature distribution for variable 
strength sources with constant spacing \ 
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Figure 3.4 Detail ofthe flat plate surface temperature distribution 
in a characteristic interval for the example of a variable 
strength - constant spacing source distribution 
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inconsistencies in the vicinity of the leading edge (as expected 
from the discussion in section 3.2.B), higher temperatures in 
regions of the heat sources with lower temperatures in the adia- 
batic regions between successive sources, and larger amplitude 
temperature variations for greater heater spacing. 

On the other hand, whereas in example 1, the temperature fluc- 
tuations increase in a streamwise direction due to the increased 
lateral heater spacing, in example 2, they decrease because of 
the reduced heat transfer. Additionally, the average plate 
temperature appears to constantly increase for example 1, but 
remains essentially constant over a major portion of the plate 
for example 2. Unlike the results of the former, where 
the temperature variation is fairly symmetrical around successive 
sources, the ununiform spatial temperature oscillations for the 
latter are thought to be caused by lateral heat transfer in the 
plate arising because of the variable source strengths. Note from 
the amplitude irregularities shown in figure 3.4 that the 
horizontal spacing between adjacent heaters is apparently not 
the only criterion for determining the amplitude of the temperature 


variations. 





CHAPTER 4 
GENERATION OF SMALL AMPLITUDE WAVE DISTURBANCES 
IN THE PLAM PLECATE BOUNDARY LAYER 
4.1 Selection of Generation Method Used 

To select the means for generating oscillatory disturbances 
in the boundary layer of the designed flat plate assembly, it is 
specified that the disturbance mechanism be capable of producing 
both "finite'' amplitude disturbances and the "infinitesimal 
Tollmien~Schlichting waves" assumed in conventional stability 
analysis. This latter capability should permit a reasonable com- 
parison between experimental results and an analytical solution 
to the problem under consideration (references (2) and (3)). 

A number of different methods have been tried in similar 
experiments. Schubauer and Skramstad (7), concerned with sub- 
sonic speeds and frequencies of 10 - 260 Hz , investigated the 
effects of sound, first by installing a loudspeaker at the 
leading edge, and then by introducing the sound into the boundary 
layer through a small hole in the flat plate; but they found that 
resonance effects and the complexity of the resulting wave patterns 
excessively complicated boundary layer perturbation with controlled 
sound intensity. They further used successfully an electro- 
magnetically vibrated ribbon stretched parallel to the surface of 


the plate and to the leading edge. This method has since been used 
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by other authors ((29) and (30), for example). Laufer and Vrebalo- 
vich (31), concerned with the laminar supersonic boundary layers 
and frequencies of 5000 - 50,000 Hz , also tried using sound waves, 
in this case generated by a cutout in a second flat plate 
(oriented parallel to the main flat plate), but the radiated 
field proved to be too wide, disturbing the boundary layer along 
the whole length of the plate. They tried an electromechanical 
device where a thin diaphram mounted flush with the surface of 
the flat plate would be made to oscillate by acting as one plate 
of a capacitor; but they found that insulation break-down in the 
capacitor at low voltages and shielding of their hot wire amplifiers 
from the high frequency voltage source radiation field made this 
method untenable as well. A third (successful) method used by them 
was essentially a high speed valve, where a narrow slit in the 
plate would be periodically opened and closed to permit an air pulse 
of any required strength and frequency to be injected into the 
boundary layer. Polymeropoulos and Gebhart (32) and later 
Dring and Gebhart (33) perturbed a natural convection boundary 
layer with a wire ribbon, but they oscillated the frame on which 
the ribbon was stretched parallel to the face of the plate rather 
than the ribbon directly. 

After reviewing all of the methods used above, the vibrating 
ribbon was chosen for this experiment. In essence the sinusoidal 


oscillation is obtained by passing an alternating current of a 
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desired frequency through the ribbon in the presence of a 
constant magnetic field. For the forces to be normal to the 
plate, the magnetic field has to be perpendicular to and in 
the same plane as the undisturbed wire ribbon. To produce 
such a field, a bipolar or horseshoe type magnet is used with 
the wire located equidistant from the two poles. An alternate 
approach might be to run a constant current in the wire and 
oscillate the magnetic field; but it was felt that this large 
alternating magnetic field would most probably induce 
voltage fluctuations in the probes used to detect the stability 
of the flow. Therefore, the first method was used. 
4.2 Compatibility of Experimentally Generated Wave 

with the Assumed Tollmein-Schlichting Wave 

At this point, the question might be raised as to whether 
such a device could closely enough model the two-dimensional 
disturbances assumed in conventional stability analysis. Gaster 
(34) mathematically shows that "a ribbon-type disturbance ex- 
cites a travelling wave system downstream of the source. This 
mode will have a wave number, & , identical to that of the free 
eigenmode which has the same frequency as the ribbon." The complex 
nature of this wave number confirms that "the mode excited by 
the vibrating ribbon is indeed one having spatial growth as 


observed in experiment." 
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Lin (35), in analyzing the work of Schubauer and Skramstad (7), 
points out the conformity of the free oscillations observed by the 
latter with the two-dimensional theory, noting that these oscilla- 
tions are known to be three-dimensional and should be expected to 
be such due to their instability at sufficiently high Reynolds num- 
bers. He shows that the two- and three-dimensional disturbances 
have nearly the same second branch point for the neutral stability 
curve and so are indistinguishable, concluding that the oscillations 
in the boundary layer become progressively more three-dimensional 
as the Reynolds number of the flow increases. Similar conclusions 
have been drawn by Klebanoff and Tidstrom (30). So, it would seem 
that any attempt to generate only two-dimensional disturbances (as 
assumed for the theoretical model in this analysis) is impossible, 
but this difficulty had been removed earlier (at least for an un- 
heated plate) by Squire (36) who proved that a two-dimensional, 
incompressible flow becomes unstable at larger Reynolds numbers 
when a disturbance within it is assumed to be three- rather than 
two-dimensional. Hence the critical Reynolds number is determined 
by the latter type of disturbance. 


4.3 Determination of Significant Parameters Affecting Vibrating 
Ribbon 


A. Explanation for analytically modelling the vibrating ribbon 
For design purposes, it is thought necessary to obtain some 

sort of analytical solution, prior to the construction of the 

perturbing mechanism, to resolve two problems. First, since such 


a device has been used successfully for forced convection studies 
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exclusively in air, where one would expect the viscous damping to 

be much smaller than that in water, prior to use in the latter 

medium an indication of the effect of viscous drag on the smooth 
oscillations of the ribbon should be examined - specifically, the 
phase and form distortion of the generated wave and change in its 
amplitude. Nondimensionalization of the governing equation with 
respect to significant parameters dispels the hope that the relative 
smallness of the coefficient of the damping term would permit this 
term to be completely excluded in an approximate solution. There- 
fore, the nonlinear viscous damping term is retained in the formula- 
tion of the problem. Secondly, it is desired to obtain an idea of 
the magnitude of the driving force needed to vibrate this ribbon. 

For these reasons the following analytical problem is posed. 

B. Analytical solution of the nonlinear, fourth order equation 
defining the motion of a "beam" string with viscous damping 
being driven by an electromagnetic force 

To identify all the physical processes taking place and to 
facilitate the formulation of the mathematical model, a statement 

of the problem to be solved is made. A wire ribbon of length 1, 

width © , and thickness K , İs stretched parallel to a flat plate 

at a close distance and constrained at each end with built-in end 
conditions (that is, slope and deflection at end points are zero). 

It is forced to vibrate electromagnetically in a viscous medium 

by passing an alternating current through the ribbon of frequency 


UD in the presence of a constant magnetic field. 
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In deriving the equation for the motion of the disturbance 
ribbon, the following assumptions are made: 

1. Motion takes place in only one plane, and in this plane, 
each particle moves in purely a transverse direction (i.e., at right 
angles to the equilibrium position of the ribbon). 

2. Deflection of the ribbon is so small at all times that 
its change in length when displaced has no effect on the tension. 
That is, the slight rotation of the cross-sections and elongations 
of the ribbon are negligible so the effects of structural damping 
(37), (38) and the spring restoring force are neglected. 

3. Since the maximum amplitude of the ribbon is assumed to 
be small, for simplicity of the solution, it is further assumed 
that there is no change in the magnetic field over the range of 
displacement. 

4. Viscous drag coefficient of the ribbon is essentially 
constant, representing almost exclusively form drag. Because 
the length of the ribbon is much greater than the width, which 
in turn is much greater than the thickness, it can be assumed that 
this ribbon acts as an infinitely long flat plate moving normal 
to a given flow. This coefficient has been determined to be a 


constant, independent of the speed of the flow (39): 


C = 2.0 


D 
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5. The ribbon is treated as a "beam" string to account for 
the bending forces. 

6. Spatial boundary conditions require that both ends be 
clamped in place: that is, the slope and displacement are assumed 
to be zero for all time. To establish temporal boundary condi- 
tions, it is stipulated that the effect of damping will be a 
maximum when the velocity of the ribbon is the greatest and zero 
when the ribbon is at rest. 

With these restrictions in mind, the equation describing the 


forced oscillation is derived. 








Sh + See Gall 
T 
Edo 2 z . 2\ 2 E st (4.3.1) 
2 
KS = IB sin Wt 
or nondimensionalizing with 
T= yt pod EE (4.3.2) 
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with spatial boundary conditions 


FREE, at r Q n 
a 
Temporal boundary conditions are to be established shortly. 
Initial significance attached to multiple-time scaling 
was soon revealed as unfounded, and so the solution to equation 
(4.3.3) was obtained by a double expansion in esos with 
the viscous damping) and E (delineating the relative importance 
of bending rigidity and ribbon tension). (See appendix C for 
further elucidation). Hence the nondimensional amplitude can be 


t= [6 ( +e (o) +] 
tt o, (e)| K (#7) p eK, (#7 q | 


Dr 


(4.3.4) 


where essentially the functions e: and Ko represent the first 
order corrections to the forced oscillation of a massless string 
due to the effects of viscous damping and bending, respectively. 
Substituting the assumed form of T into equation (4.3.3), an 
ordered set of equations in powers of E results for zeroth order 


in € . Imposing the spatial boundary conditions and applying 
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the intricate matching conditions in regions near the ribbon's 
end points (where it is confirmed that the uncorrected solution 
must satisfy only the homogeneous spatial boundary conditions - 
that is, temporal boundary conditions need not be considered 
until the effects of damping are introduced), the zeroth order 


solution in both e and € is found to be 


Be (E,7) > £ | tan = SiN AF + Cosas = Isin r (4.3.5) 


Note again the Lo is expected to be the major contribution to 
the solution for the ribbon's displacement. 

With go thus determined, the constant coefficients in 
the expression for (6) ( z T ) can be found by satisfying 
the matching conditions and by realizing that the viscous 
damping must be zero at the time when the local velocity of the 


ribbon is zero. Therefore, 


gi (% 7) = = [K (g) ses z + K, (5) Cos 37 


+KP cos OT + K cos Ie tal CoD 
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where 


> 
, (tang SIN & + Cos a _ 1) 


3 | | 

K _ i tan S Sin Sat + Cos Sx f =| 

2 45 (tan SE u oe, | 
oa 


K Su K p Sin De F + COS Sa id 
3 MA A AR + cos Ar = 1 


VASE: 
K e ¡[E EN + cos Taf = 
2 
Again applying matching conditions, it can be shown that for 


first order matching in € , Y (€) = le and so 


K, (8,7) = Ltang|- tan$ sınaf - Cosaf| sin (4.3.7) 


Extending this procedure further, K, and additional correction 
terms can be obtained, but, for purposes of this analysis where 
even the first order corrections are shown to be small, it is 
felt that higher order approximations are unnecessary. 

In an attempt, now, to obtain some idea of the effect of 
the viscous damping and bending on the ribbon's movement, 
typical values for the problem are used to evaluate the solution 


(see table 4.1), 
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Material: spring phosphor bronze ribbon (94.75% Cv 3% Sn, 23% P) 





K =ı01. 

K = .0035 in 

Ke = ,006 in 

oO = ‚0625 in 

m = .37 x 10° slugs/ft 

Ce = 17.05 Sie Fr. 

E = 16 x 10° psi 

Í = 9 in 

L = IE = 6.51 x 10°" in’ 
oe 


1.94 ae 


(z> 
I 


20 


Table 4.1 
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The significant parameters of the problem can be determined: 


e = 2.06 x/0°” 

&”= 1.765 

@ = eo 

% is postulated to be of unit order. 


Hence, the time dependent amplitude of the ribbon at its midpoint 


is 


t= | . 1530 SsınT (.00073 COS, OOO Cosa 


8 | o oo we) 


-|. 00340 sin T + O (pre) | 


From this result, it is seen that the effect of the damping 

is to cause a damped ribbon to lag an undamped one (i.e., ee 0 
solution). Additionally, since the ribbon is required to complete 
one cycle of its oscillation in the same period as that of the 
driving force, it is verified that the damping will tend to speed 
up the ribbon in some regions (immediately after the maximum 
displacement is reached) forcing it to go slower than the undamped 
one in others (near regions of zero displacement). These two 
expected effects are qualitatively displayed in figures 4.1 and 
4.2; quantitatively, the effects are much smaller than those shown 


as substantiated by the above sample calculation. 
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Effect of damping and bending on the displacement 
of a vibrating ribbon. 
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Figure 4.2 Effect of damping and bending on the velocity of 
a vibrating ribbon. 
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C. Simplified solution to the same problem 

Since it has been shown in the preceding section that the 
effect on the wave form of a vibrating ribbon due to damping and 
bending is small, for subsequent analysis, higher order correc- 
tions are neglected and the displacement is approximated by 


the zeroth order solution. That is, 


E E (8,7) 


= Y 
x? 


(4.3.5) 
[tang sin x? + Cos xf = ] SAA 


or expressing all terms dimensionally, 


yes sin (o f)2 + sin( 213) 


| un wt 
7 sin( RM > (4.3.8) 
K Mm 
- Sin (w an l 
which corresponds (as it should) to the solution obtained by 
Hussain and Reynolds (29) where the bending and damping force 
terms were not even considered. This simplification of the 
solution verifies that the central portion of the ribbon 
essentially behaves as a massless and completely flexible vibra- 


ting string, satisfying the equation, 


5 . í 
N: = K = + ] BSinwt 
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with conditions, 


y (0, ) = y (1,t) = O 


D. Attainment of the proper mode and amplitude of oscillation 


To this point, no attempt has been made to examine the 
character of the resulting solution for ribbon displacement (i.e., 
points of resonance, number of nodes, effect of ribbon tension, 
etc.) or to answer the second equestion posed in section 4.2 
regarding the magnitude of the force required to drive the ribbon. 

Examination of equation (4.3.8) reveals that resonance 


would occur when 


or ae N” 
n 1 | m 


where w, is designated the natural’ frequency of oscillation. 
Since it is desired that the ribbon oscillate in a single 

loop standing wave as visually verified in use in other experiments 

(7, 29), parameters affecting the number of maxima/minima must be 

adjusted until this configuration is attained. By maximizing the 

equation (4.3.8) (and sine that 7/1 < 1 ), it can be shown 


that multiple loops form only if 
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Ole t= ern F 412 3; 
(Then = 0 case is just the desired single loop.) Apparently, 
then, multiple loops should not occur until after having passed 
through the second resonant condition. After that point, for 
2tr< WE Lede three loops (or four nodal points) character- 
ize the motion, for 47r< I I< Cir five loops should 
theoretically be observed, and so on. Figure 4.3 shows the first 
two modes of oscillation. 

With the formula for 4 ( Ze ), the force required to drive 
the ribbon at a certain frequency for a desired maximum amplitude 
can easily be determined, keeping in mind that w [2 | Zi 
for a single-loop oscillation. This amplitude, occurring at the 


ribbon midpoint can be found from equation (4.3.8) to be: 


oT | - Cos w me 
m (t) j ea —— SUN wt 
COS CO a 2 


s a tan (af 2) tan[e| = 2) sin wt 


For example, assuming 








les 10 Aen ps. = .37x 10 Sie 


10 lb 
= qin 


rad 


m 
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Figure 4.3 Significance of the parameter ol f% in determining 


the mode of oscillation of the disturbance ribbon 
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the magnetic flux intensity, B , can be determined, 


=m n pu. 





Note for this example, oj] = .9064 < 27 and so only a single 
loop should be expected. 

Just as the maximum and minimum allowable temperatures for 
the existing testing configuration were computed in section 
3.2A, limitations are imposed on the range of frequencies required. 
According to Wazzan, Okamura, and Smith (reference (2), table 2 
and figure 30), the nondimensional maximum frequencies of 
unstable disturbances for the previously specified temperature 


interval (0 £ je = 1, € 22.5°F, T = 60°F) are 


aD 


23 
a = /D 
Unheated plate TOS max 
= o = _4 
where U) = WY 
May sS 


2_ 
Uco 
Again using as the limiting tunnel speeds, un 40 fps and 


US 5 fps, 
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5980 Hz 


maximum frequency 


9,4 Hz 


minimum frequency 


4.4 Implementation of the Method for Vibrating the Ribbon: 


Since other experiments which used a vibrating ribbon to 
perturb the flow were conducted in a wind tunnel, very little 
difficulty was encountered in putting a given amount of tension 
on the ribbon ~- using rubber bands and weights. In this design 
where problems with sealing against water leakage from the tunnel 
were encountered, special provisions had to be made to counter 
these problems while still permitting adjustment of the ribbon 
tension. Rather than using suspended weights, the design shown 
in figure 4.4 was used - a glorified spring-loaded fish scale! 
Note that while tightening of the bolt in the plexiglas frame 
to compress the packing (which sealed around the ribbon) will 
seal the hole, it will also tend to give erroneous readings on 
the calibrated tensioning device if adjustment is attempted while 
the packing is compressed. 

This is not felt to be too great a drawback, since tension- 
ing could be made prior to sealing tightly. Once the tension is 
set, it should not need to be readjusted. In a similar line of 
thought, the edges of the plate would have to be sufficiently 
rounded at the points where the ribbon would pass So as not to 


excessively bend the latter. 
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Figure 4.4 Disturbance ribbon tensioning device 
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To eliminate possible electrical current paths other than 
through the wire ribbon (a potential shock hazard), the 3/4" 
outer cylindrical part of the device was threaded into a teflon 
insert which had been press-fit into the aluminum supporting block. 
This cylinder, in turn, was encased in 5/8" x 3/32" clear tygon 


tt 


tubing, and tensioning was accomplished by a teflon "nut" at the 
top of the cylinder. A previously noted precaution, specifically 
anodizing the plate to create a dielectic surface, eliminated any 
path from the wire to the plate (see appendix A-2). Since the 
packing device is imbedded in the nonconducting plexiglas frame, 
no current leakage should exist to the water tunnel. 

The inability to easily vary the strength of the magnetic 
field was sacrificed for the elimination of the power source 
and temperature control devices required for the proper operation 
of electro-magnets (see (29) for example). Hence, two sixty pound 
pull, ALNICO, permanent magnets were chosen. Since it was realized 
that the magnetic flux density would drop rapidly away from the 
magnet face, a recess was milled in the 3/8" thick flat plate to 
place the magnets to within 1/16' of the working side of the plate. 

With other design considerations in mind, it was determined 
that the closest that the ribbon could be placed to the leading 
edge was 2.5". These ribbons were made from spring phosphor 
bronze due to that material's high fatigue resistance, tensile 


strength, and electrical conductivity. They were initially .1' 
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wide and .0035" thick, and stretched across the width of the plate, 
9", It is anticipated that perhaps the width, length of vibrating 
span, and distance from the plate might need to be adjusted to 
produce a measureable disturbance, although Schubauer and 

Skramstad (7) indicated that for operation in a wind tunnel, none 
of these factors significantly affected the performance. Naturally, 
the further into the boundary layer the ribbon remained, the less 
flow- and self-induced flutter should be expected. Although the 
amplitude of these disturbances could be analytically calculated 
from the previously discussed model, it is thought that the ampli- 
tude constituting an "infinitesimal" disturbance could only be 


determined by experimental measurement. 
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APPENDIX A 


ANODIC COATING ON THE ALUMINUM BASE PLATE 


A.1 Formation of Anodic Coating 


The anodizing process used was the ALCOA Alumilite coating 226, 
standardized as 2 mil thickness. The properties of the dense, ad- 
herent anodic oxide coating to be discussed in section A.2 are a 
function of the electrolyte solution composition (a 12% sulfuric - 
1% oxalic acid solution), bath temperature, composition of the 
alloy, the current density, and the elapsed time. As pointed out 
by Wernick and Pinner (20), the actual mechanism of anodic oxida- 
tion is very complex and controversial, and dependent on the types 
of electrolytes used,to determine the composition of the coat- 
ings. In simple terms, in this electrochemical process an electric 
current is passed through an aqueous electrolyte solution in which 
the aluminum piece to be anodized acts as the anode and the lead 
walls of the container act as the cathode. The negatively charged 
anions consist in part of oxygen ions which react chemically 
with the aluminum to form the anhydrous aluminum oxide in accordance 


with the formula: 


2 Al de 5 O —ALO, + Energy 


Since the reaction products are sparinply soluble in the electro- 


lyte, the film growth on the surface of the anode is accompanied 
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by its dissolution as well. It should be noted that the oxalic 
acid is added to the electrolytic solution to decrease the solu- 
bility of the coating serving to increase its hardness and 
limiting thickness. This dissolution forms pores in the coating 
which are wide enough to allow continued access of the current 
to the metal (specifically, the pore diameter is on the order of 
S O celle inches 23) ieetielincandefiim stones 
retarded by the increased electrical resistance caused by the 


thickening film, continues until the rate of growth and dissolu- 


tion are equal, at which point the film thickness remains constant. 


The coating will also contain small percentages of sulfate and 
oxalate. In addition to the primary reaction, a number of 
secondary reactions also occur at the anode such as the reduction 


of oxalic acid: 


(CooH), + 0 —— 2C0, + H,O 


and the electrolysis of water (40). Gas analysis at the anode 
has indicated the presence of carbon dioxide, carbon monoxide, 
oxygen, hydrogen, and nitrogen with proportions depending upon 
me current density (20). 

After plugging desired holes with rubber or plastic plugs 


to restrict the area to be anodized, the piece was then put in 
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the electrolyte, ensuring positive electrical contact, and the 


) 
current turned on raising slowly 2-5 volts per minute to the 
prescribed current density (in this case 36 EN Due to 
the discharge of energy in this reaction as indicated above, the 
electrolyte had to be agitated and refrigerated (8-10°C) since 
insufficient agitation would cause locally higher electrolyte 


temperatures, which, in turn, could cause burning. The thinner 


the sheet, the more agitation was required to prevent burning. 


A.2 Effect of Coating on Plate Performance 


An assessment of the various properties of the anodized 
coating and how they would affect the performance of the plate 
was necessary. A Summary of that study is presented below. 
Hardness and abrasion resistance 

An assessment of the wearing quality as measured by con- 
ventional indentation tests are not indicative of the material's 
hardness, since such test points load the coating, forcing it 
into the much softer Al substrate producing too deep an impression 
for reliable Fes litem For example, repeated sharp blows from 
a hammer can dent and slightly chip the coating, owing to the 
thinness of the coating and the softness of the basis metal. 
However, we have already specified that the plate le not be 
subjected to localized loading, so we turn to the results of 


two abrasion tests - resistance to wear by rubbing and resistance 
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to corrosion from a grit blast. In such tests the resistance to 
abrasion of the hard anodize coat was found to be about twice 
that of the regular coating (169 gm/mil of coating to 89 gm/mil 
for std coating, 13 times that of untreated Al, and comparable 
to case-hardened steel and chromium plate. For rubbing wear, 
6061 T6 hard coated Al performed about like case hardened steel 
and much better than carbon steel or type 304 stainless steel. 
Resistance to corrosion 

Even though the pores formed in the hard anodizing process 
are not Sealed, since that would reduce the abrasion resistance, 
the coating provides excellent resistance to atmospheric and 
saltwater corrosion. Compared with conventional coatings, which' 
are usually expected to withstand a few hundred hours of salt 
Spray exposure without failure, 6061 alloy with a 2 mil hard 
coat withstood 10 months exposure as per ASTM B117-49 with no 
sign of failure (15). Experiments on the effect of rain erosion 
on the leading edges of high speed airfoils, conducted by Gillig 
show that failures depend primarily on the alloy which is hard 
coated. The best results were obtained with 24S Alclad, 61S, 
and 24S, in that order (14). 
Fatigue and tensile strength 

The anodizing process causes a slight decrease in the ten- 
sile strength (for 61S T6 with a 2 mil coating a reduction from 


47,700 to 47,000 psi) believed due to the stress concentrations 
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at the microcracksin the coating. In general, since anodic 
coatings in excess of .5 mil tend to reduce the fatigue strength 
of wrought Al, hard anodic coatings are generally not used when 
the fatigue strength is a factor; but for this application, it is 
not a factor. The slight loss in ductility is not considered 
important. 
Adherence 

Due to the nature of the process, the coating becomes an 
integral part of the piece. The only question of adherence would 
be at the leading edge as discussed above. Although severe bend- 
ing of the piece should not be experienced, tests show that a 
right angle bend over a 3/4 inch diameter shows that a 2 mil 
coating spalls off on the compression side of the bend and checks 
finely on the tension side, but the piece can be bent 180 degrees 
on a liberal radius and straightened out again without spalling 
off of coating. The adhesion of the oxide coating is normally 
much better than that of electrodeposits. 
Electrical properties 

Since the wire ribbon of the perturbing mechanism could 
come into contact with the plate, it would be advantageous to have 
the latter insulated. For the nonconducting, unsealed anodic 
coating (2 mil), the breakdown voltage is 950 volts. Although 
sealing in boiling water of the coating increases the breakdown 


voltage for 2 mil coating to 1,200 volts, it is not considered 
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necessary to obtain the additional resistance at the expense of 
the abrasive resistance (although the plate would still be well 
within the limits required for this application in either the 
sealed or unsealed configuration). 
Effect on surface finish 

Due to the excellent throwing power of the process, uniform 
coatings may be applied to irregular-shaped objects but varying 
slightly for a closed end hole (i.e., the thickness of the coating 
in the unplugged blind bolt holes in the plate could not be 
expected to be uniform). In general the surface roughness is slightly 
increased by hard anodizing. For example, for 7075 T6 with a 
coating thickness of 1.2 mil, an initial surface smoothness of 5-7 
microinches (RMS) will have a smoothness of 10-20 microinches 
after the treatment. If a finer finish were desired, the plate 
could have been made oversized and then lapped or honed down to 
a smooth surface on the order of 2-3 microinches. Since the plate 
was machined, then lapped prior to treatment, it is felt that 
the surface smoothness is well within the limits at which it 
will contribute to turbulence generation. 
Heat resistance and thermal expansion 

In general, due to the differences of coefficients of thermal 
expansion between the anodic coating and the aluminum substrate, 


checking of the film occurred when the temperature exceeded 200°F. 
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The coating itself acts as an insulator - a disadvantage when 
trying to conduct heat to the surface of the plate to be con- 
vected into the boundary layer - having a thermal conductivity 
of 7.17 Btu/ (hr) ft” °F/in. (41). The coating does, however, 
have a good resistance to flame impingement or short term high 
temperature exposures. When heat transfer is by radiation 
emissivities of 90-95% are to be expected. Due to the highly 
refractive nature of the coating, the interior will melt before 


the coating collapses. 





APPENDIX B 


ANALYTICAL SOLUTION OF HEAT CONDUCTION IN A COMPOSITE 
MATERIAL WITH HEAT GENERATION IN ONE REGION 


on me ee ee ee eee PO 


Heat is generated by a finite series of rectangular heat 
sources in one region of a two-region, composite rectangular slab 
which is insulated on three sides but which transfers heat by 
convection to a laminar boundary layer on the remaining side. 

It is required that the temperature distribution on the plate 
surface adjacent to the flow be found so that the influence of 
spacing and strength of the sources can be evaluated. 

B.1 Determination of Source (Neumann) Function 

Following the notation of the mathematical model shown in 
figure B.1, the equation to be solved with accompanying boundary 


conditions can be formulated: 


K (EN au): Ko Y Ni = - S(x- 5 (y- 4) (B.1.1) 


(where the point source discontinuities are in region da ), with 


boundary conditions 


aloe Delay SE Gy = Se (a,y)70 











(B132) 








7 ie) >. Aad e 
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Note in the specification of boundary conditions that the Neumann 
function(s) is (are) derived with homogeneous boundary conditions 
(see (42) for elucidation). 

Bulavin and Kashcheev (43) have solved the transient heat 
conduction problem, for multilayer symmetrical bodies with heat 
generation within the solid, by the method of separation of variables, 
and by construction of orthogonal expansion of functions over multi- 
layer repions for the case of perfect thermal contact between 
the surfaces. Similarly, for the steady problem here, it is 
assumed that N. X ' J ) can be expanded in an orthogonal 


eigenfunction expansion, 


N. (y 5x4) -)) n N Gey’) Yin (J) cos I 


mzo n=f 
i 
on PR (B.1.3) 
Sian =! vn Z O 


where a are eigenfunctions satisfying the following differential 


y” Yin(y) a O, Yin (y) =a) (B.1.4) 


Now considering the orthogonality property of the eigenfunctions 


i. over the entire rectangle, i.e., over the ] layers, where 
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“| ato Ye (y) 7 (constant) Sn (BS) 


with ' being the steady, discontinuous weighting function 
(thermal conductivity) for layer i which is unchanged for either 
the assumption of perfect thermal contact or linear contact 
resistance at the interface or for boundary conditions of the first, 


second, or third kinds (44). 


A b go 
I N. cos T ao ne) I NY cos"z dydx 


= ag Ky | Ya dy 


| Jj 





1 


(B.1.6) 


jet y? 
where M= K: (B.1.7) 


(Bee) 


NY) = 2 | DA N Yi. cos I dy dx 
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Applying integral transform (8.1.8) to equation (B.1.1), 


Zo e) f [Ne EN a Tax V Ade 


a yo 2 
ON. EN » TMX 
En J | (2 Ne a 0) cos tm Y. dudx 


(B'1 0) 


“ a 
| 9 Na C o j. dy = aNa Cos MX 
x a Ox A 





A 


a , fà 
mm Si Tm) (mm TMX 
ò Ó 


or after applying boundary conditions from equation (B.1.2) 





A 2 a 
3 Na mx = 
cose = Le) cos TMY d (B.1.10) 
| Ox? a NA s 
° ° 
Similarly 
A : a a 
O Na ces OS dx =- = ea In Cos TMX dx (B.1.11) 
D A A < a 
© O 
also 


„g 








ü 








b b 
Na EN i er Yan 
s C: : dy af dy Nad e 
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or applying equations (B.1.2) and (B.1.4) 


= > Na Y dy = ze (x, D) 











(B.1.12) 
EP 0) = To - Ji hA ay 
eo 
BEN, e de 5 a (D SDK (x, o) 
aC al Je A (Bol: 13) 


- N. eae i d e (o) - ee JY on Noy 
Substituting equations i 1.10) through e 1.13) into equation 


(5 1.0), 


ae: a: 
= EA + On lU | CN en dy dx 








Jfk NY, Cos z dydx! Ail eG > (50) 





PEN Lalo + OFEN 





an (o daai: = — | 
= K.N.&o) = )| Cos —— Arf T coa mx Yan? 


or upon simplification with equations (B.1.2) and (B.1.8) 
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Tx 


Cos 
N yl zn (B.1.14) 
=D 
with equations (B.1.3) and (B.1.14) thus established, the 
problem is essentially solved; all that remains is to determine 


the expressions for Yin € q and On . Recalling equation 


(B.1.4) where 
Le . L £ 
Ç Yin (4) a Ya: O 
it can be shown that the solution takes the form 
Yin (4) = A., Pin (Y) + Sin Un (y) 


where On q and Vin a are two linearly independent 
solutions, which in rectangular coordinates take the form 

Pin (y) = sin Rey 

Win (y) ES a 


Therefore, 


O j Pe Sin <j b. COS ig (B.1.15) 





115 


Since Yin Su ) must satisfy boundary conditions in equation 


(B.1.2), equation (B.1.15) can be substituted into these boundary 


conditions to yield 
An COS y - Bn Sin nb =© 


Tat Bs = 0 


(B.1.16) 
WER: 
A. Cos pac + B.. Sin Bae = O 
DS a Cot Ba b 
A g a Aa (B.1.17) 
Den” >. A an cot En D 
and 
lt Ka ctac = QO (B.1.18) 
ë K. On 


Hence, Rn 


( y ) is determined in terms of one constant 
(here chosen as An which will cancel as the solution progresses, 


and On are the eigenvalues determined by equation (B.1.18). 
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From equations (B.1.15) and B.1.17) 


” CoS Bn(y-D 
N ü Din Sin Ab 


Yan? -Aan Ka cos eny +9 


RS 


(B: 1.19) 


Now, to determine = from equations (B.1.7) and (B.1.19), with 


os = a e + = 








2 + 6. 
and 
w G Nor A Cm . 
J Yen g | > 5 * Cn 
then | 
i Yin 2 
Mz) KY, (y) dy 
a N 
PRA, 


Hablan lbsiWenc + Kac sineb 
2 sun"@,b seek š ae + ie Sin E | 


(B.1.20) 
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Combining equations (B.1.3), (B.1.14), (B.1.18), (B.1.19) 
and (B.1.20), the source functions are completely determined, and 


may be written as 





fT COS al -b) Sa a a(y-B) Sin! E C Cos stm E 
Na ~ ak, = ren JL b sià enc + Fee sin'enb | 


mao n=1 
and (B.1,21) 
on cO I = 
| ` M MY 
N,= 4 Im COS - b) cos (y+c)sın b sın cos eos a. 
Eee ak. 


m0 n=l E ) +e Ife Sin a= + Ke sine, b | 


where $ = 1 m 4 O 


B.2 Determination of Temperature Distribution 

With the source functions now determined, attention is directed 
toward finding the temperature distribution on the surface of the 
rectangular slab resulting from a series of rectangular heating 
elements of strength ss thickness E , and width A , lo- 


cated at positions ( Xe E Ji ) in region A . To represent these 
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sources, a combination of step functions is used so that the 


internal source distribution is 


° H| y- (4:5)! - Aly u (y:- Sy (280) 
S, (x,y) =0 


where it is assumed that Ss is uniform within each element, 
and that the dimensions on all elements are identical. Using 
the procedure outlined by Duff and Naylor (42), with the kernel 
or Neumann source function now defined by equation (B.1.21), 
the temperature can be determined to within an additive con- 


stant with any desired internal source distribution and speci- 


fied surface heat flux as follows: 


O-T-Tae= | Nis dA + N. pi 44 (B.2.2) 





17 


where a 


local temperature 


average temperature over the region - £ [Taa 
| A 


internal source in region | 


D I on the external boundaries of nu ] 
J on 
subject to the restriction that in a steady state configuration 


the heat generated internally must be transmitted through the 


boundaries. 


2 

fs aA = -|v N;dA = = oN dl = ~ { p dÍ (P.2. 3) 

À À Pen + 

If this condition is not satisfied, then no solution can exist. 
Recalling the discussion of the model to be used in this 

analysis as presented in the introduction to this appendix, the 


following boundary conditions are stipulated: 


u 9- -— 
i po (B.2.4) 


P(x% y) O on all other boundaries 





EN 
where for want of a better choice of a constant le K 2 wo Al 
Y 


in accordance with the similarity, Blasius-type solution for a 


thermal boundary layer. 
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Incorporating eqns. (B.1.21), (B.2.1), and (B.2.4) into 
equation (B.2.2), the temperature distribution in the plate 


would be 


CO e I | = 
= E Sa) Cos = ( E un = “ee - b) T En 
a L — x b sin 2 E 0 Sin On 


L = 
| a. 

Be? ) Ces Pe TE) Sin = sin eis Cos | H dy’ 
bs eine] J T 
n=1 Ñ | Ñ Ke € O 

a TMY 
3 WN BAC 
aS Fee tet ones aa 
i — LEN + On bsin Ge rc enb] 


o E i Sa as ln un u cospaly; - bco, (y-ysin bak 











al os Ë, Sin ab cs @ (De fai sim (B.2.5) 
KS SB. (y+c) N cos P| 19) Be ee 


a prodigious looking equation to say the least! 
Furthermore, to insure the validity of a solution, the con- 


dition specified by equation (B.2.3) must be satisfied. Therefore, 
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J En ‘| pe (1-94 
A 2 
or ) SDE AA (B.2.6) 


After attempting to determíne the eigenvalues with the charac- 


teristic constants, it becomes apparent that the eigenvalues are 


rn 
approximately En TE <= "Uli Sp ecicicallys 


feeb = .25", C= .002", K = 117 Btu/hr ft F, and K = .597 


Btu/hr ft F, 


@, = 3.999830 = Bods 


T2) 


e = 7,999673 
2 25 


N 


It can be seen that the real eigenvalues will be less than the 
approximate ones (a difference that will increase as N—>0 ), 
but owing to the relative rapidity with which the solutions are 
expected to converge, this difference is not expected to signi- 
ficantly affect the solution. 

Y 


\ 
Therefore, specifying that On = —— + the source functions 


of equation (B.1.21) are drastically simplified, becoming, 
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m cos Tmax Gog mx 


ES 4 ) Y cos 4 aa 
= abK, mn U) + TO | 








N. © | (p. 2.7) 
h 
where un Í meo, Paa 

re Eu n#o 

¿a 2 M#O 4 n=O 


os 30 


Physically, this implies that the lateral heat transfer in region C 
is insignificant when compared with that in region A . Consequently, 
the slab should be treated as a single component material and the 
source function, or kernel N, be used with all sources and 


boundary conditions in the integral equation (B.2.2). 


_ 4 cos TAY cos MEY sin Ht ma 
Ka. D I 4) Sa P i mey = 
L | n=1 a 


co 
Ten X am A 
oe Cos ` cos EN e ya mm dx’ 


a L, e 
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een Ta 

cO 
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This equation can be simplified somewhat by evaluating some of 


the series. From Jolley (45) 


oo ee C Ta (b- | 
en ee lm 


and 





co 


Er = a b“ 
EY |Z 
Making the appropriate substitutions, the temperature distribu- 


tion is finally obtained. 


Bi g A cos MAY cos MY! sin THE 
le gp 
a7 ES mr 


E, y E A E (B.2.8) 
m 


mS Eo “s sb s sin ME 
L. (=) w (GS 


m=1 n=l 





> A 
RO b E 2 cosh UH (6- Y, Tim \ H Tx 
|a: a Fad T ca sinh SE osm [cas s— dx 
m=1 
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B.3 Numerical Examples of Temperature Distributions for a One 


eee 


Material Moc _ Model for Specified } Heat Source Di Distributions 











As an example of the solution to equation (3.2.9), the required 
heat flux was assumed to be provided by twelve sources. An initial 
indication of the strengths and positions of these sources is ob- 
tained from equation (3.2.14), derived in section 3.2C, for the 
following configurations: 

l. Constant strength sources with variable spacing 

2: Variable strength sources with constant spacing 
Pertinent parametric values applicable to both are found in 
table 3.1. Source strengths and positions in the plate for con- 
figurations (1) and (2) above are found in tables 3.2 and 3.3, 
respectively. 

In the evaluation of equation (B.2.8) for a specific example, 
the Fresnel integral is encountered, i.e., an integral of the 


form 
im 


ss cos Ta = DG 2 2 
. dx = (EE) cos u du 


which must be approximated as follows (see Jolley (45)): 


. 1)" (Tm) = for small TTM 
Er, Bo 


° 
— em | = |. 3:95 Pu for large IT m 
| YZm' (2itm) (201) | | 
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In practice the first formulation is used for m=3 and the 
second for all higher values of M . 

To obtain the desired accuracy in the temperature distrubu- 
tion, k 6 of equation (B.2.8) is summed on the computer so 
that the individual infinite series are converged to an accuracy 
within3.01. The FORTRAN V program used is shown in figure B.2. 

Surface plate temperatures resulting from the example con- 
figurations (1) and (2) are shown in figures 3.2 and 3.3, respectively, 
with characteristic details of the latter appearing in figure 3.4. 

A discussion of the ramifications of these examples and figures 


is presented in section 4.2D. 
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APPENDIX C 


ASYMPTOTIC SOLUTION TO VIBRATING RIBBON 


C.l Formulation of the Problem 

A wire ribbon of length 1, width @ , and thickness K , is 
Stretched parallel to a flat plate at a close distance and 
constrained at each end with built-in end conditions (i.e., the 
Slope and deflection at the end points are zero). It is forced 
to vibrate electromagnetically in a viscous medium by passing an 
alternating electric current through the ribbon of frequency u) 
in the presence of a constant magnetic field. To determine the 
effect of small viscous damping and bending on the motion of the 
ribbon, the above model is formulated and solved with the desired 
effects being just perturbations on the main forced oscillation. 

Consistent with the notation in the mathematical model 


shown in figure C.l, the defining equation is 


(Gals) 

LY, = mh + ehe (ape E Kp = IBsin at 
— / Bu, Ee tot ñ j 

bending inertial viscous u. forcing 

force force drag force function 

Nondimensionalizing with 

Zr 12 La £ = — 
= l 
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Wire Ribbon 


Mathematical model used to examine the electromagnetically oscillated wire ribbon 


Figure C.I 
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equation C.1.1 becomes 




















st. of 2 af] 
€ T sint C.1.2 
oF >F? + Ir In? Í In Ir CA ( ) 
where 
é = Elo Y = 181° 
r Ki 


T S T 
with boundary conditions; 
spatially: ¿ir at on 
temporally: the requirement that the viscous damping be a 
maximum when the ribbon velocity is a maximum 
and zero when the ribbon is locally at rest. 
Note that the resulting equation is a fourth order, nonlinear, 
partial differential equation. Examination of this equation 
suggests two possible ways of reaching an approximate solution: 
1) a regular perturbation in Ce > solving the fourth order 
equation exactly, or 2) a boundary layer type solution, with 
interior and exterior expansions matched in terms of an inter- 
mediate variable. The interior solution would be expanded in 


a regular perturbation expansion in N ; 
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The first approach leads to such a complex zeroth order solution 
(no damping) with the end effects incorporated within it, that 
the higher order corrections (where the effect of the damping is 
considered) prove to be too cumbersome for our purposes here. 
Specifically, since the behavior of the wire ribbon near its 
center - away from the end points - is the quantity of interest, 


the second method suggested above is used. 


C.2 Asymptotic Solution for Interior Region 


With the viscous damping, it is expected that perhaps there 
might be more than one time scale associated with the problem; 
that is, one associated with the damping ( Tr = eT ) and one 
with the oscillation of the wire (= Í + e a) ) where W, 
is a frequency shift caused by the damping. 

Following Cole (46), it is assumed that for the interior 
portion of the ribbon, the dependent variable of the defining 
equation can be expanded as 


t < m ty OLT (Fone S 8) 


(2) 
* Z. 
O 1 +... 
where WV (ce) y (¢)... are some functions of € to be determined 
l > ì ) 
later. By substituting equation (C.2.1) into equation (C.1.2), 


the ordered set of equations is powers of € is obtained. 
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a 2 
Ë Sp uw = SE 5 sin (ca. 


(Note that the time used in the forcing function is that asso- 


ciated with the oscillation of the wire.) 





St 
= h. pye 1 
2 e L = 
ch aa, =e |. = j (C.2.3) 
š a l 
OF OT ar f oo 


Further, for small damping, define 


h.= Po ee + OQ: Ms 7) = e (E MT)... 
hy K ET un) + oK G m) tee. (C.2.4) 


so by substituting equation (C.2.4) into equation (C.2.2) (for 


the zeroth order solution in € ) an ordered set of equations 


in £ is obtained. 


É Ba o = A “ee 


ml 
o 


Er 
g: 2 








133 


By considering the zeroth through the second order equations, 
utilizing the fact (determined later by matching the interior 
and exterior solutions - see equations (C.3.7) and (C.3.15), 
matching conditions (1) and (5)) that the zeroth order solution 
in E must satisfy homogeneous boundary conditions at £ SD. Es 
and requiring that the solution be uniformly valid in pr. it 
can be found that ¡lo through Oz are not functions of ar 
and that U), = 0. This means that “T = Y ; in other words, 
there is only one time associated with the problem. In 
retrospect, this seems reasonable for,since the wire ribbon 
is driven externally, its frequency should be unaffected by 
the damping. 

Therefore, equations (C.2.1) and (C.2.4) can be redefined 


and combined to obtain (to first order in ( ) 


{= h, (?,7) t Y, (€ ) h, (£, 7) +... 
= | En) + BB (FT) vera +- Y (€)| Ko (¥,7) +0 Kent. 


a (Cyne) 


Equations (C.2.5) for zeroth order in € can now be rewritten as 
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O: A s. 











9 E am? 
e: $ en + q 2 br == do (G.2 7) 
Dp? ar For 
2 _ See 25 e aa >e, PRI 
C ore Bre ` 


It follows, then, that the zeroth order solution in 7 and € is 


Bo (E,7) = sin TLA, sino F + 6, cost - x 


+ tos C, Sina E i D, cos | (G.258) 


C.3 Asymptotic Solution for the Exterior Regions 
and Matching with That for the Interior Region 


To determine the exterior expansions, it is assumed that the 


bending and the tension terms are of the same order of magni- 
tude. Stretching the Ç variable in the boundary layer area 
near Ç = 0 so that 

~ r 

(A 


the governing equation (C.1.2) becomes, 








O is at it, we 


an = 1: 
DE Sog OTa oT 





vr 
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By equating the coefficients of the bending and tension terms, 


the exponent,  , is found to be n a= 1/2. 


Therefore, the 
boundary layer coordinate near r O is 


Sd 


2 = (C.3.1) 
Further, expanding 


E 15€) L en OREN HORE neil 


the governing equation near f = 0 becomes 


E pos Wu Ep = 
HS” ie ot E +098) 


af lok 
KA nr ISE 
OT 


ar 
To lowest order, two choices are available for the form 


1. Emo 
































or 


as E — 00 


(09.29 
— 


Assuming (1) is valid (until contradicted or confirmed) 


an ri 
at — af 6 
mA 


T: (0.359) 
>š oF 





» 


By satisfying the boundary conditions at f = 0, and assuming 
= 
that the solution of the form € will not match with the 
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nd 


interior solution as f O 4 
paaa T ME 
f, = A Lt ME ] (0.3.4) 


Proceeding formally, an intermediate variable is chosen 


suitable for matching near £ = 0; i.e., it is desired that 
PS 
de while Bu as. G6 


€ — O 
where fr is fixed and for { 
) 


sie 


Defining + = 
1 


ne 
Fog oh — 


IS 


as € —O (C.3.5) 


In preparation for matching the interior and exterior 
solution, equations (C.2.6) and (C.3.4), the former is expanded 


in a Taylor series expansion about Ë = 0. 
f = Le (0,7) m Ñ ey (O+. : We f Le. (0,7) + eg (o) T) 
De 0. | + E LO” (o T) ER On e ol Reh 
gr Y ol Ka + e K (o T)+...| r v,(€) [LK Cr) 
e CK (ott. rt 
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where primes indicate differentiation with respect to spatial 
coordinate, f 

Matching the interior and exterior expansions in terms of 
the intermediate variable Í near £ = 0 | 

] 

lim E 
ne [ee + CB (0,7) de le: (o T) + e gr(o, 7) 
F, tixed 
taJ Ht LK EDA OKC +++] 


+ En AO LK.) +e K Co, T)+. | +... = (€) Alber, 


Ye 


> 
-{+e ano 1 (Ask, D | = O Cc.3.6) 
Matching is possible only if: 
1. Py (0,7) + 89,00 += 0 
2. Me (€) = YE 
(6.347) 
3, A BolT) te RODE.. 


n AE [KO eK N.a] 


Note that since Mule) Y € , the assumption that FE > is 


justified. Note further that the term of 0(1) in E cannot be 
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matched except by a term from h , implying that 


1 
OIR (Gere) 


Since WE) u E—06, the zeroth order solution in O ín the 
Ye 
first order solution in € can be found (see equation (C.2.3)). 
Following the same procedure at the other end of the wire, 


the stretched coordinate is defined as 
et 
E = (C.3.9) 
yé 
and the expanded deflection i becomes 


f= (ef (gt re... (310) 


Putting this in the governing equation (C.1.2), and applying 


boundary conditions, 
+ + i 2 
[ = A mlite BER: 


Now defining an intermediate variable so that near € = 1 


“tP 
as i. | i u GZ 2 C€ — Ó 





r = EOD 


g (€) 


SO (C.3.12) 
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o~ (€) — O G 
where ee is constant and as 7 a j Ta 


Therefore, 


a a) —— | 
(G; 3.15) 


T(E) 
o = Verte 7 78 





as E —=0 and T(E)— 0 wi th fr fixed, 
As before, matching of the interior and exterior solutions 


in terms of the intermediate variable is accomplished near Ç = 
lim / 

+6 Mr re [tet [e | K. EKA Dto 

EYE e (ALKOT) + eK Cr) +...) Fe 


PNE TE ie - 
- YE AFA Fe tl-e Jade (C.3.14) 
Matching is possible only if: 
5: Bol T)+ ee COT) t= 0 
en An: es (17) = °G Ga... (C.3.15) 


7. No Zee eK, n 
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As noted previously when showing that only one time scale 
is associated with the problem, the zeroth order solution in € 
must satisfy the homogeneous end conditions (matching conditions 
(1) and (5)). 

With the main contribution to the solution for the vibrat- 
ing ribbon, Bo , now specified by equation (C.2.8), first 
order corrections due to the influences of viscous damping 
(represented by E ) and bending (represented by Ko ) can 
be evaluated by the solution of equation (C.2.7), and equations 
(C.3.8), (C.2.3), and (C.2.4), respectively, with matching and 
boundary conditions determining the constants. 

C.4 Utilization of Matching and Boundary 

Conditions to Determine the Complete Solution 

Utilizing the information derived in the preceding two 
sections, the first order corrections to an undamped, completely 
flexible oscillating ribbon due to viscous damping and bending 
can be found. 


Bo - the basic solution 





Using matching conditions (1) and (5) for equation (C.2.8) 
and realizing that each term in an asymptotic series should in- 


dividually satisfy the boundary conditions, 
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Qo (5,7) = (A, sink + ID, Cosa F = e) Sin T (0.208) 
+ (C, Sinx E + D, cos ar) COS 


so that 
P(O T) = . Sin T +D, om. (ENS 
Lo Klo) = | A, sins +B, Cose ls T+ Cosina J: Dies STE O 


To be valid for all Y 


_ ` [| 3 Goce) E Y A 
A, MS a O = 


O 
| 


9 |< 


C. © (Ga) 
D, O 

so 

= ca T) = Stan sinaf + tosaf- ı|sinr (C.4.3) 
nr ci 


With the form of Eo now determined, equation (C.2.7) can 


be solved to yield EN ( r A s To Simplify the Solution 
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the nonlinear damping force, being an even function in T 


is expanded in a Fourier cosine series 


de. |de. a 
aa — > b, Cos NT 
Ma 








Or (C.4.4) 
2 
= tanz Sin xe + COS x € z | COS [Cos T| 
where = 
L 2 X 
b = RE Ba ss: cost = (| | “cos Cos ny dr 
Aa z ) 
H: 2 
- | Be Neo ar 
T 
de 
Or 


2 
“alle 
> ES tang Sun AF + COS AF Pi en 


sn (2-N)% 4 Sin(2e)Z] n#2 (6.4.5 
SO) 2(2+n) 
Therefore, 
b = 0 for n an even integer 
= 
b, = ° Lo tan sina + cos«f -1] 
25 
2 57 (0.4.6) 
Le 
b; * "z£ 
- 
> 105 


Using equations (C.4.4) and (C.4.6), equation (E 


may be 
written as 
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2 
ee - b, (Cost + Ll cos 37 
= 1 r~a 1 © o ` (C.4.7) 
== Cos Š + 155605 IT a) 


Applying the matching conditions of equations (C.3.7) and 


(C.3.15), which essentially require that the ribbon satisfy 


homogeneous conditions at its end points, equation (C.4.7) yields 


pi (pI = Br (tanz sin af + cesu- i) Cost’ +, (tan $ sin 3uf 


nn 1) Cos ST rt Side E ai cos Jak - 1) cos dr 


+ aoe ( tan Te to Sin Ing + COS (aE = 1) Cos [Tr -.. .| 


or 


gT) = “KW Cost + EY a 
+ K, (Í)tos 97 K (cos 71 te (C.4.8) 


where 


es; 3 
( = [tant sima? + COS af -1) 
K ... Kt) (AE + o 


5 Y tans ; F 


K.G) =- (E ) = SınSx + ccS S< - j 


e75 + COS AF - | 


í / tan lt sin Tx} + Cos Tay - +) 


as tans sino $ + Crs FI 


oe 
et 
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K bending correction to basic solution 


with f ( a , T ) now evaluated to zeroth order in € , 


a first order correction due to bending is sought. Since the 


form of V, ( € ) is determined in equation (C.3.8), equations 
(C.2.3) and (C.2.4) become (after revision necessitated by 


the presence of only a single time scale), 

















O 2 = 
~ - DKo 4 y? 28 =0 
97? MN 
i: p Sk, Mea > K. JE u (C.4.9) 








From zeroth order equation, Ko can be found to take the 


following form 


e (A, sina i O, cosa) Sin + (C,sine E +D, Cosa $) Cost (C.4.10) 


From matching conditions (3) and (4) 


er one. =- Ko(0,7)+ eK, (O T)+,,. 


so that, using equations (C.4.2) and (C.4.10), for zeroth order 


in e 


= (0,7) = A, Salz S, sin T o D ces T on 
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Similarly, using matching conditions (6) and (7), 


| IE Zee Kur re Kant 


so for zeroth order De 


ay) = Cw sa = B xsinx )sinT = (A sing = 


O, Cos) SsinT + C, sina + D, Cosa) Ce ass s js) 
To be valid for all time, 


A, =~ fA, tang) = - 2 tang) 


ee - AA, = “x tang) 
C= 0 
DAO 


(C.4.11) 


Consequently, with equations (C.4.10) and (C.4.11) 


KF, a "Stan Etang) sinag + Cos | ET (6:94.12) 
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